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ABSTRACT 


The combinatorial geometry allows efficient Monte Carlo 
radiation transport calculations of detailed three-dimen¬ 
sional geometries. Because regions and media are formed 
by combination of basic bodies such as boxes, spheres, 
cylinders and others, the input required of the user is both 
relatively simple and easily modified. The M0RSE code is 
a multigroup neutron and gamma ray transport Monte 
Carlo code that may solve either neutron, gamma ray, or 
coupled neutron-gamma ray problems in either the forward 
or adjoint mode. M0RSE has a wide variety of available 
input options, including splitting, Russian roulette, expo¬ 
nential transform, energy biasing, importance regions, 
albedo surfaces, and the scoring options available in the 
SAMB0 analysis package. 

This document details the incorporation of a new version 
of combinatorial geometry into the M0RSE code and is 
meant as a user's manual. All new subroutines, including 
the combinatorial geometry, are documented in detail, in¬ 
cluding flow charts. Input instructions for all M0RSE 
options and layouts of all M0RSE commons are given for 
user reference. The PICTURE code has been modified to 
work with the combinatorial geometry, allowing two-dimen¬ 
sional printer pictures of combinatorial geometry models. 
Two M0RSE sample problems and a PICTURE sample 
problem are discussed. Input and output from these prob¬ 
lems and the complete code are available from the 
Radiation Shielding Information Center. 
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I. INTRODUCTION \ 

The M0RSE Monte Carlo code was first distributed^ in 
early 1970 with a revised version^’ released in the fall of 1970. 

There have been several additional features^’added in the past year. 
This report describes the addition of a version of the MAGI developed 
combinatorial geometry' ' (CG) to the M0RSE code and also serves to 
document a version which is operational on the UNIVAC 1108. The 
definition of variables in labelled commons, layout of blank common 
and input data have been collected for easy reference. A description 
of the PICTURE program for the combinatorial geometry package and 
descriptions of routines required to solve two sample problems are 
included. \ 

The replacement of the 05R geometry module with the CG 
module required significant modification of the CG package. These 
modifications were required so that the CG piodule could be added to 
M0RSE rather than adding MORSE to the CG package. Only those 
M0RSE routines that have an interface with the geometry package were 
changed and all options in M0RSE are available with the CG module. 

To retain all M0RSE options it was^ necessary to write several new 
routines as part of the CG module. ( 

V Features of M0RSE which are dependent on the geometry 

module include the ability to (1) deter mine, both the, media and an im¬ 
portance region given a set of coordinates, (2) to track a particle through 
the system including internal voids and surrounded by an external void, 

(3) determine the number of mean free paths between any two sets of co¬ 
ordinates (one may be in either an internal or an external void), (4) 
determine the vector normal to a surface at any point, and (5) permit 
particle reflection at a boundary without "getting lost." Table I lists 
the M0RSE module routines which were modified to permit an interface 
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TABLE I. 


M0RSE Routines, that Depend on Which 
Geometry Module is Used 


EUCLID (3) 
G0MST (2) 
INPUT 1 
MS0UR (1) 


with the; CG module. The numbers after the routines Relate the features 
mentioned above with the function of the routines. The routines without 
numbers are concerned with the input of geometry information. Table n 

TABLE II. 


Geometry Module Routines 


05R Spherical Geometry 

Combinatorial Geometry 

GE0M (2,3) 

G0MFLP (5) 

J0MIN 

^ L00KZ (1) 

N0RML (4) 

GTV0L 

V 

GENI 

GTVLIN 

GG (2,3) 

G1 \ (2,3) 

ALBERT 

PR 

G0MFLP (5) ; 

J0MIN 

L00KZ (1) 

N0RML (4) ' 

GTV0L 
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lists the geometry module routines for the 05R spherical and the com¬ 
binatorial geometry packages. 

Only minor changes were made to INPUT and MS0UR and no 
additional documentation is provided. INPUT was split into two routines, 
INPUT1 and INPUT2, to allow a more efficient overlay. Numerous 
changes were required for subroutine EUCLID and G0MST and new de¬ 
scriptions and flow charts of these routines are given. A complete 
description of the combinatorial geometry package is included and 
sample problems are provided. 


3 










II. 


COMBINATORIAL GEOMETRY 


2.1 Body Types 

Combinatorial geometry (CG) describes general three dimen¬ 
sional material configurations by considering unions, differences 
intersections of simple bodies such as spheres, boxes, cylinders, etc. 

In effect, the geometric description subdivides the problem space into 
* 

unique zones. Each zone is the result of combining one or more of 
the following geometric bodies. 

1. Rectangular Parallelepiped (RPP) 

2. Box (An RPP randomly oriented in space) 

3. Sphere 

4. Right Circular Cylinder 

5. Right Elliptic Cylinder 

6. Truncated Right Angle Cone 

7. Ellipsoid 

8. Right Angle Wedge 

9. Arbitrary Convex Polyhedron of 4, 5, or 6 sides. 

Body types 2-9 may be arbitrarily oriented with respect to 
the x, y, z coordinate axes used to determine the space. Body 1, a 
special body described below, must have sides which are parallel to 
the coordinate axes. 

The basic technique for the description of the geometry con¬ 
sists of defining the location and shape of the various zones in terms of 


To avoid confusion between importance regions and combinatorial 
geometry regions, we depart from previous combinatorial geo¬ 
metry descriptions and use the term zone to indicate a combinatorial 
geometry region which is designated by the variable IR. The term 
region is reserved for an importance region. Thus the zone index 
is IR. 
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the intersections and unions of the geometric bodies. A special oper¬ 
ator notation involving the symbols (+)> (-), and (0R) is used to describe 
the intersections and unions. These symbols are used by the program 
to construct information relating material descriptions to the body de¬ 
finitions. 

If a body appears in a zone description with a (+) operator, it 
means that the zone being described is wholly contained in the body. If 
a body appears in a zone description with a (-) operator, it means that 
the zone being described is wholly outside the body. If the body appears 
with an (0R) operator, it means that the zone being described includes 
all points in the body. In some instances, a zone may be described in 
terms of subzones lumped together by (0R) statements. When (0R) 
operators are used there are always two or more of them, and they re¬ 
fer to all body numbers following them, either (»•) or (-). 

Techniques for describing a particular geometry are best il¬ 
lustrated by examples. Consider an object composed of a sphere and a 
cylinder as shown in Fig. 1. To describe the object, one takes a spher¬ 
ical body (2) penetrated by a cylindrical body (3) (see Fig. 1). If the 
materials in the sphere and cylinder are the same, then they can be 
considered as one zone, say zone I (Fig. lc). The description of zone I 
would be 

I = 0R + 2 0R + 3 

This means that a point is in zone I if it is either inside body 2 or inside 
body 3. 
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If different materials are used in the sphere and cylinder, 
then the sphere with a cylindrical hole in it would be given a different 
zone number (say J) from that of the cylinder (K). 

The description of zone J would be (Fig. Id): 

J = + 2 - 3 

This means that points in zone J are all those points inside body 2 
which are not inside body 3. 

The description of zone K is simply (Fig. le): 

K = + 3 

That is, all points in zone K lie inside body 3. 

Combinations of more than two bodies and similar zone de¬ 
scriptions could contain a long string of (+), (-), and (0R) operators. 

It is important however to remember that every spatial point in the 
geometry must be located in one and only one zone. 

As a more complicated example of the use of the (0R) operator, 
consider the system shown in Fig. 2 consisting of the shaded zone A and 
the unshaded zone B. These zones can be described by the two B0X's, 
bodies 1 and 3, and the RCC, body 2. The zone description would be 

A = + 1 + 2 
and 

B = 0R +3-1 0R +3-2 

Notice that the 0R operator refers to all following body numbers until 
the next 0R operator is reached. 







Fig. 2. Use of 0R operators. 

The geometry must be specified by establishing two tables. 
The first table describes the type and location of the set of bodies used 
in the geometrical description. The second table identifies the physical 
zones in terms of these bodies. The input routine processes these 
tables to put the data in the form required for ray tracing. Because the 
ray tracing routines cannot track across the outermost body, all of the 
zones must be within a surrounding external void so that all escaping 
particles are absorbed. Also no point may be in more than one zone. 

The information required to specify each type of body is as 

follows: 

a. Rectangular Parallelepiped (RPP) 

Specify the minimum and maximum values 
of the x, y, and z coordinates which bound 
the parallelepiped. 











z 



Fig. 3. Rectangular Parallelepiped (RPP). 
b) Sphere (SPH) 

Specify the vertex V at the center and the scalar, 
R, denoting the radius. 



Fig. 4. Sphere (SPH). 

c) Right Circular Cylinder (RCC) 

Specify the vertex V at the center of one base, 
a height vector, H, expressed in terms of its 
x, y, and z components, and a scalar, R, 
denoting the radius. 


„Ymax 
—Zmax 
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Fig. 5. Right Circular Cylinder (RCC). 

d) Right Elliptical Cylinder (REC) 

Specify coordinates of the center of the base 
ellipse, a height vector, and two vectors in 
the plane of the base defining the major and 
minor axes. Presently this body is not 
implemented. 




Fig. 6. Right Elliptical Cylinder (REC). 








e) Truncated Right Angle Cone (TRC) 

Specify a vertex V at the center of the 
lower base, the height vector, H, expressed 
in terms of its x, y, z components, and two 
scalars, R^ and R 2 , denoting the radii of the 
lower and upper bases. 



Fig. 7. Truncated Right Angle Cone (TRC). 
f) Ellipsoid (ELL) 


Specify two vertices, VjandVg, denoting 
the coordinates of the foci and a scalar, R, 
denoting the length of the major axis. 



Fig. 8. Ellipsoid (ELL). 







Specify the vertex V at one of the corners 
by giving its (x, y, z) coordinates. Specify 
a set of three mutually perpendicular 
vectors, ai , with aj and a 2 describing the 
two legs of the right triangle of the wedge. 
That is, the x, y, and z components of the 
height, width, and length vectors are 
given. 


h 

v 

Fig. 9. Right Angle Wedge (WED). 

Box (B0X) 

Specify the vertex V at one of the corners by 
giving its (x, y, z) coordinates. Specify a set 
of three mutually perpendicular vectors, a^ 
representing the height, width, and length of 
the box, respectively. That is, the x, y, 
and z components of the height, width, and 
length vectors are given. 
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Vx, Vy, Vz 


Fig. 10. Box (B0X). 

i) Arbitrary Polyhedron (ARB) 

Assign an index (1 to 8) to each vertex. For 
each vertex, give the x, y, z coordinates. 
Each of the six faces are then described by a 
four-digit number giving the indices of the 
four vertex points in that face. For each face 
these indices must be entered in either clock¬ 
wise or counterclockwise order. 








2.2 


Subroutines 


Descriptions of the combinatorial geometry routines are given 
with logical flow charts depicting their functions. Two M0F&E routines, 
G0MST and L00KZ required significant modifications for use with the 
CG package and descriptions of these routines are also provided. 

2. 2.1 Subroutine G1(S, MA, FPD, L0CREG, NUMB0D, IR0R, 

mi, ifi2) - 1 - 1 - 

G1 is the control routine for the combinatorial geometry. On 
one call, it calculates the distance travelled in the present zone, and 
the number IR of the next zone to be entered. Essentially GG is called 
for each body adjacent to the present zone, calculating RIN and R0UT, 
the distances to entry and exit of the body along the trajectory. The next 
zone to be entered is determined by again calling GG to calculate RIN 
and R0UT for each body adjacent to the next possible zone. These next 
possible zones are determined by examining a list cf all the previous 
zones entered on crossing this body. RIN and R0UT are checked against 
the input zone descriptions to determine the correct zone. If it is not 
found in the list of previous zones, all other zones are examined in a 
similar fashion, and when the correct zone is found, it is added to the 
list of previous zones for that body. If the new zone is different from 
the old, G1 returns; otherwise G1 continues tracking until a different 
zone is encountered. One change added to the M0RSE version of G1 is 
that if the distance to the next boundary is greater than the distance to 
scattering, G1 returns without determining the next zone past the bound¬ 
ary, setting the flag MARKG in common 0RGI. 

Called from: G0MST, EUCLID, MESH 

Subroutines called: GG 

Commons required: PAREM, G0ML0C, DBG, 0RGI 
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Variables required: ! 

XB(3) - starting coordinates of present trajectory, 

WB(3) - direction cosines of trajectory, ^ 

IR - present zone, 

DIST - present distance from XB(3), 

DISTO distance frbm XB(3) to next scattering point, 

NASC - less ihan zero if this is d new trajectory, 

'j KL00P - trajectory index, 

PINF - machine infinity, 

MA, FPD, L0CREG, NUMB0D, IR0R, ' 

. IR1, IR2 - locations in blank common used for variable 
dimensioning. 

Variables changed: , 

KL00P - trajectory index incremented if this is a 
new trajectory, 

NASC - next Ipody intersected by trajectory, 

LSURF - I surface of body NASC crossed at next inter¬ 

section (negative if leaving and positive if 
^ entering NAS0), 1 

>DIST ' - distance from XB (3) to next intersection or 

collision site, ^ * 

MARK - set to 1 if distance to collision/ (DISTO) is 
\ less than distance to next intersection (other¬ 
wise 0), 1 

S - distance travelled on this call to Gl, 

IRPRIM - zone to be entered on boundary crossing, 

MA\lNEXT) - new zone added to the list of next possible 
zones for body NB0, 

MA(INEX) - location in MA of next item in the next possible 
zone list for body NB0 (these lists leap-frog 
through the end of the MA array). 
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\ 


Significant internal variables: \ 

NB0 - absolute value is body being considered 

while a negative or positive sign indicate^ 

, that zone IR or IRP is outside or inside the 

body respectively, 

R0UT - : distance to exit of body NB0 calculated by 
GG, 

RIN - distance to entry of body NB0 calculated by 

I GG , 

LRI - surface of body NB0 entered by trajectory, 

LR0 - surface of body NB0 trajectory exits, 

! IRP - zone being considered as next zone. 

Limitations: 

I01JT - output logical unit is set to 6 in data state¬ 

ment. 


\ 
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2.2.2 Subroutine L00KZ (X, Y, Z, MA, FPD, L0CREG, 

NUhffi0D, IR0R, NS0RT - - 

The purpose of this routine is to return the combinatorial geo¬ 
metry zone of point (X, Y, Z) so that tracking can be initialized. The 
coding has been borrowed from the second half of subroutine G1 and 
adapted to determine the zone of a source particle. For efficiency 
L00KZ builds a li3t of possible source zones to search on future calls. 

If the region is not found on this list, all other zones are examined and 
upon determining the new source zone, it too is added to the list. Notice 
that the starting direction cosines (.8, .6, 0.0) are assumed in L00KZ, 
but may be changed elsewhere. 

Routines called: GG 

Commons required: 0RGI, PAREM, G0ML0C, and DBG. 

Variables required: 

X, Y, Z - coordinates for which a zone is desired, 

MA, FPD, L0CREG, NUMB0D, IR0R, NS0R 

locations in blank common used for variable 
dimensioning. 

Varaibles changed: 

KL00P - trajectory index is incremented, 

NMED - common 0RGI variable set to correct zone 
number, 

NS0R(INEXT) - new zone added to list of possible source 
zones. 

Significant internal variables: 

WB(3) - set to . 8, . 6, 0.0 so that L00KZ need not 

be called with a direction 


CG version. 
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Subroutine L00KZ (X, Y, Z, MA, FPD, L0CREG, NUMB0D, IR0R, 
NS0R) 


/ Is position identical N_ Y 

V to last point? / 


Increment KL00P and initialize 
XB(3), WB(3) and DIST. 


Do ove r 20 possi b le source zones IRP . 


(Is IRP zero? )-Yes- 


*0 


Loop over all bodies adjacient to IRP. 


- (Is NB0 negative? 


Is R0UT < 0 or \ ( Is RIN<DIST 

R0UT < DIST -no—- no— and 

or RIN > DIST?/ V DIST<R0UT? 




I 


¥ 


NMED « JRPRIMl 


fRETURl 









m 


_ T 

Save index for adding a zone 
to the source list 


Loop over all zones IRP. 


Loop over all bodies defining zone IRP. 










2.2.3 Subroutine GG (L0CAT, MA, FPD) 


GG is the workhorse of the combinatorial geometry, comput¬ 
ing distances to intersections for all body types. It is called from G1 
or L00KZ to compute distance to entry and distance to exit to a body 
whose location in the MA array is flagged by the argument L0CAT. Each 
time a distance to entry, RIN, and a distance to exit, R0UT, are calcu¬ 
lated for a body they are stored in the MA array together with LRI and 
LR0, the indices of the entry and exit surfaces of the body. Also stored 
at this time is KL00P, the particles trajectory index. On a subsequent 
call to GG for that body, KL00P is checked against the earlier value, 
now L00P. If they are the same, the old values of RIN, R0UT, LRI, 
and LR0 are retrieved so that GG can return immediately. 


If it is necessary to compute a new trajectory a different area 
of coding is entered for each body type to calculate RIN, R0UT, LRI 
and LR0. 


Called from: G1 and L00KZ 

Subroutines called: none 
Commons required: PAREM, GE0M2 
Variables required: 


LOCAT 

KL00P 

PINF 

MA, FPD 


starting locati n in the MA array of integer 
data for the appropriate body, 

trajectory index, 

machine infinity (stored in RIN and R0UT 
when the trajectory misses the body), 

locations in blank common required for 
variable dimensioning. 
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Variables changed: 


RIN 

distance to entrance, 

R0UT 

distance to exit, 

LRI 

- surface of entrance 

LR0 

surface of exit. 

Significant internal variables: 

L 

same as L0CAT, starting location in MA 
array of body data for the appropriate body, 

K 

starting location in FPD array of floating 
point data for the appropriate body. 
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Subroutine GG (L0CAT, MA, FPD) 


Retrieve K, L00P, and 
IT YPE from MA. 


_j_ j Set LRI, LR0, RIN 

Has this trajectory previously \ y ( and R0UT from MA 

v been calculated for this body ? ' '— A TT,r * T '- 


and FPD arrays. 


-non Write error i 



Store KL00P in MA array! 
for this body. _| 


( Transfer on IT YPE to compute RIN, R0UT, LRI, and LR07) 
(AREl) (SPH)(6 cC?) (REC) (TRC) (El£) (B0X) (RTI 


Store LRI and LR0 in MA array. 


Set RIN or R0UT to DIST when 
one is within roundoff of DIST J 


iR0UT=R0UT*l.OOOOl| 


|R0UT =-infinity! 


Write 
’’Type not 
implemented’ 


I Store RIN and R0UT 














2.2.4 Subroutine GENI (MA, FPD, L0CREG, NUMB0D, IR0R, 

BRIE, MRCZ, MMIZ, MMCfe, KR1, KR2, il, I0UT71N) 

Although subroutine J0MIN is used to read in all of the com¬ 
binatorial geometry data except the region volumes, that data is put into 
the proper storage location in blank common by GENI. This is accom¬ 
plished by having J0MIN write the body and zone data on a mass storage 
unit, so that GENI can retrieve that data. GENI assigns the data to 
blank common in the area set aside by J0MEN. For an arbitrary poly¬ 
hedron body, GENI calls ALBERT to handle a portion of the data. In 
addition, GENI outputs the input information concerning the geometry 
data. GENI computes certain geometry dependent data needed during 
the random walk. 


Called from: 

J0MIN 

Subroutines called: 

ALBERT 

Commons required: 

Variables required: 

Blank, G0ML0C, PAREM 

11 

- input unit, 

I0UT 

output unit, 

IN 

bulk storage logical unit. 

Variables input: 

Geometry data for zones and bodies is retrieved 
from mass storage unit IN. 

Variables changed: 

MA(I), FPD(I), L0CREG(I), IR0R(I), MRCZ(I), 
MMCZ(I), KR2(I)--see Fig. 13, layout of com¬ 
binatorial geometry storage in blank common. 
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Subroutine GENI (MA, FPD, L0CREG, NUMB0D, IR0R, MRIZ, MRCZ 
MMIZ, MMCZ, KR1, KR2, IN) 











2. 2. 5 Subroutine ALBERT (F, IERR) 

ALBERT is called by subroutine GENI to process the arbi¬ 
trary polyhedron (ARB) body data before storage in the FPD array. The 
ARB body data as read by GENI consists of the coordinates of all corners 
(eight for a six-sided figure), followed by a packed decimal number for 
each side indicating which corner points make up that side. ALBERT 
processes this data, returning a unit normal vector and a minimum 
distance to origin for each plane containing a side of the ARB. An ARB 
can have up to six sides. These unit vectors and distances then replace 
the original ARB data in the FPD array. The number of sides and a 
distance that is characteristic of the ARB's minimum dimension are 
also stored in the FPD array. This minimum distance is later usi... for 
round-off tests. 


Called from: GENI 

Subroutines called: none 
Commons required: none 
Variables required: 


F 


FPD array read from GENI 


Variables changed: 

F - FPD array changed to contain unit normals 

and distances to origin for each plane. Num¬ 
ber of sides and a minimum distance are 
also stored. 


Significant internal variables: 

X(8,3) - coordinates of cornei points, 

IS(6,4) - indicates corners contained in each of the 

possible six surfaces, 

V(3,4) - edge vectors for a given surface (an inner 

cross product then gives the normal vector). 


29 








Subroutine ALBERT (F, IERR) 


[s tore coordinates of corner points in X(8,3) and zero F array.| 


ICount the sides' and store ^ corner indices in IX(6, 4),| 
|Calculate minim um dist ance between points!] 


("Begin a calculation for each sideTj 


[C ompute two adjacent edge vectors. 1 

_u_ 

Compute normal vectors as the cross products 
_ of two adjacent edge vectors. _ 


Compute niinimum distance to orgin as dot product of I 
normal vector and coordinates of the corner.' | 

I 

j Count thenumber of points on either side of the plane. | 


(A re all points on th e far side of the plane?)- 


IChai.ge sign of the normalizing factor, j 


Normalize normal vectors to unity. Also 
I normalize the distance by the seme factor. 

^- I: 

[lsjjhis the las t side?") - no- 

1 yes; 

_ .{ _ : _ 

Store normal vectors, distances, minimum side lengtn 
and numocr of sides in FPD array. 


1 
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® Subroutine PR (Combinatorial Geometr y) 

Subroutine PR is called from various locations in the combi¬ 
natorial geometry package (GENI, Gl, L00KZ and N0RML) whenever 
intermediate or debugging output is required. The amount of geometry 
data which is output depends on the value of the argument, ‘if this argu¬ 
ment is 1 or 8, all of the geometry data in blank common is printed, 
otherwise only selected variables are output. By comparing the argu¬ 
ment value given in the output with the source listing, the geometry data 
at a, given time in the ^execution can be determined. The call to PR is 
initiated by setting the IDBG variable to a nonzero valued Because of 
the large amount of output generated, this option should not be used dur¬ 
ing a normal execution. 

Called from: GENI, Gl, L00KZ, N0RML 

Subroutine called: none 

Commons required: blank, PAREM, G0ML0C, DBG 
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0: K, XBG), TOO), IR 

IR, IRPRIM, NASC, LSURF, LEU, LR0, KL00P, 
L00P, ITYPE, N, NUM, L0CAT, ISAVE, 

INEXT, IRP, INEX, LDATA, RIN, R0UT, SMIN, 



0: G0ML0C Array 
MA Array 
FPD Array 

L0CREG, NUMB0D, IR0R, MRIZ, MRCZ, MMIZ, 
MMCZ, KRl, KR2, and NS0R Arrays 
Region Volume Array 










2.2.7 Subroutine N0RML (MA, FPD, L0CREG, NUMB0D )* 

The purpose of subroutine N0RML is to return a unit vector 
to a combinatorial geometry body NASC at point XB + D ‘ * V?B which 
must be on the surface LSURF of body NASC. This unit vector is useful 
either for albedo scattering or boundary crossing flux estimates. The 
sign of the unit vector is chosen so that • ii is negative meaning that 
the unit vector will point against the particle direction. 

Called from: ALBD0 

Commons required: N0RMAL, G0ML0C, PAREM 

Variables required: 

NASC - body number particle is on, 

LSURF - surface of body NASC r 

XB(3) ) 

WB(3) > - trajectory parameters, 

DIST ) 

MA, FPD, L0CREG, NUMB0D 

locations in blank common used for variable 
dimensioning. 

Variables changed: 

UN(3) - direction cosines of unit vector stored in 

common N0RMAL. 

Significant internal variables: 

XP(3) - set to XB + DIST * WB, 

X(3) - usually used to relate the intersection point 

to a body centered coordinate system, 

H(3) - usually used as a body orientation, 

L - body location in MA array, 

K - body location in FPD array. 


CG version. 

























2.2.8 Subroutine r~VUN (FPD, MRIZ, VN0R, IV0PT, NIR, 

NtZ, 11, 10; ----- 

This routine is called by J0MIN to read in or to calculate the 
volume of each region in the geometry. The four options available are 
(1) to set each volume to 1., (2) to calculate the volumes for concentric 
spheres, (3) to calculate the volumes for slabs (not coded at present), 
or (4) to read in the volumes for each region from cards. The volumes 
are stored in blank common and are only used by the track length and 
collision density estimators. 


Called from: J0MIN 


Subroutines called: none 
Commons required: none 
Variables required: 


FPD 

MRIZ 

VN0R 
IV0PT 
NIR 
NIZ 
II, 10 


array containing zone description data, 

array which relates M0RSE region to input 
zones, 

array containing the M0RSE volumes, 

options for determining volumes, 

number of regions, 

number of zones, 

input and output logical units. 


Variables changed: 


Variables in blank common starting at KV0L. 
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Subroutine GTVLIN (FPD, MRIZ, VN0R, IV0PT, NIR, NIZ, II, 10) 
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2. 2. 9 Subroutine J0MIN (NADD, II, 10) * 

The input of the geometry data is controlled by the J0MIN 
subroutine, which performs the following tasks: 


• Reads all geometry input data except the region volumes. 

• Writes the body and zone data on a mass storage unit 
(I0UT-16). 

• Determines the length of all geometry arrays. 

• Calculates the beginning location in blank common of 
geometry arrays. 

• Initializes geometry arrays. 

• Calls the GTVLIN subroutine which returns region 
volumes. 

Since combinatorial geometry input data is dynamically allocated to con¬ 
serve storage area, it is stored temporarily on a mass storage device. 
This allows the core storage requirements to be determined. Hence, 
much of the coding in J0MIN is similar to GENI, which reads the data 
on the mass storage device and puts it into blank common. 


Called from: 
Subroutines called: 
Commons required: 


INPtfTl 

GENI, GTVLIN 

B^ank, G0ML0C, PAREM 


Variables required: ^ 

4 

All variables* in G0ML0C, II, 10 

4 

Variables input: 


IV0PT, IDBp, NAZ, MRIZ(I), MMIZfl) 

A 

Variables changed: A 


All variable^ in G0ML0C, IV0PT, IDBG, MRIZ(I), MMIZ(I). 
A 

* A 

CG version. 
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Important internal variables: 


I0UT 

NAZT 


mass storage unit, 

total number of adjacent zones summed over 
all zones. 
















2. 2. 10 Subroutine GTV0L (MXREG, GN0R) 

This routine is called by ENDRUN to calculate the reciprocal 
of the volume of each region. It is geometry independent, since it uses 
information determined in GTVLIN. 


Called from: ENDRUN 

Subroutines called: none 
Commons required: G0ML0C, blank 
Variables required: 


KV0L - index from common G0ML0C, 

MXREG - number of regions in the geometry. 

Variables changed: 


GN0R 


reciprocal of the volume ior e ion. 
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2. 2 - 11 Subroutine G0MST (TSIG, MARK) * 

Any boundary crossing between the present and next collision 
sites are determined by calling the combinatorial geometry routine Gl. 
Before the Gl call, combinatorial geometry variables in common PAREM 
are initialized, and after the call NUTR0N variables are updated. 

WATERS is summed in blank common as a track length flux estimator. 
MARK is set to -1 for an external void and -2 for an internal void. 

Called from: NXTC0L 

Subroutines called: Gl 

Commons required: NUTR0N, AP0LL0, PAREM, 0RGI, G0ML0C 
Variables required: 

X, Y, Z, U, V, W, NMED see common NUTR0N for definitions. 


BLZNT 

value of IR from last track or from L00KZ, 

XB(3) 

starting coordinates of present trajectory, 

WB(3) 

trajectory direction cosines, 

MARK 

flag to indicate type of trajectory, 

ETATH 

distance to be travelled in cm if the flight 
remains in the same media, 

DIST 

present distance from XB(3). 

Variables changed: 

X, Y, Z 

endpoints of flight, 

DISTO 

distance from XB(3) to next collision site, 

ETAUSD 

mean free paths travelled on this call to 
G0MST, 

ETATH 

cm travelled on this call to G0MST, 


CG version. 
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MARK 


NMED and 
NREG 

BLZNT 


flag indicating type of termination of flight, 

0 - normal boundary crossing, 

1 - flight within one medium, 

-1 - particle escaped, 

-2 - particle entered an interior void, 

medium and region of end point, 
combinatorial geometry region of end point. 
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Subroutine G0MST (TSIG, MARK) 










2.2.12 Subroutine EUCLID (MRK, XI, Yl, Z1, X2, Y2, Z2, P1P2, 
IG, ARG, NT, MEDIUM, IBLZ, NREGN) * 

\ \ 

Subroutine EUCLID is used to determine the number of mean 
free paths between two points. However, if the argument NT is non¬ 
zero only the mean free 'paths to the first media boundary will be calcu¬ 
lated. The routine works much like NXTC0L and G0MST. After ini- 
tializihg the combinatorial geometry variables, the particle is tracked 
from point one to point two. Total cross sections are determined by 
calling NSIGTA and distances to next interfaces are calculated by call¬ 
ing Gl. On an initial call to EUCLIE), MRK should be 1. If NT is non¬ 
zero, MRK should be 0 on successive calls for the same trajectory, so 
that new trajectory parameters are not initialized. , This will be handled 
automatically if the calling routine does not change MRK. 


Called from: \ 

GETETA, RELC0L 

Subroutines called: 

Gl, 

NSIGTA 

Commons required: 

Variables required: 

PAREM, 0RGI 

MRK 

- 

set to 1 for initial call of a trajectory; 

XI, Yl, Z1 

- 

coordinates of starting point, 

X2, Y2, Z2 

- 

cooi^jdinates of end point, 

' P1P2 

- 

distance between starting and end points. 

IG ' 

- 

energy group index, 

NT 

- 

0 for total mean free paths, 4 0 for mean 
free paths between intersection points, 

MEDIUM 

- 

NMED of point one, \ 

NREGN 

- 

NREG of point one, i 

IBLZN 

- 

IR of point one. 


CG version. 
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Variables changed: 
MRK 


XI, Y1,Z1 

ARG 

NT 

MEDIUM 


1 for a flight reaching the end point, 

0 for a flight crossing a medium boundary 
(NT 4 0 only), 

-1 for a flight escaping the system, 

-2 for a flight encountering an internal void 
(NT 4 0 only), 

returns boundary intersection point if 
NT 4 0, 

negative of number of mean free paths, 

if NT 4 0 on input, will return as -1 if an 
escape occurs, 

medium number of end point. 


Significant internal variables: 

ETA - distance remaining to point two in cm, 

ETAUSD - distance travelled on last call of Gl. 






Subroutine EUCLID 











2.2.13 Subroutine ENDRUN 


This subroutine is called at the end of a run. The weight of 
the particles at collision and the track length of the particles has been 
summed over all collisions within the region for all particles in the run. 
Thus, a normalization by the number of particles in the run is required. 

In addition, the volume of the region and the total cross section for each 
region and energy group is used to determine the average fluence over 
the regions for the collision density estimate. For the track length esti¬ 
mate the same coding is used with the total cross section set to 1. Note 
that a region cannot contain more than one medium . Since storage al¬ 
locations are flexibly dimensioned, NEX in the SAMB0 input must provide 
for MXREG + 2 arrays in Blank Common to be used by ENDRUN. A 
check is made and if this condition is not met, no output from ENDRUN 
is obtained. 

There are many manipulations in ENDRUN with the energy 
indexes. These manipulations are required in order to treat the six 
major options in M0RSE. Appendix A gives an example of the six prob¬ 
lems. 

Called from: NRUN 

Subroutines called: GTV0L, NSIGTA 
Commons required: blank, PDET, USER 
Variables required: 

NMED (REGI0N) - must be set in data statement, 

GN0R - reciprocal of the volume for each 

region. 

Many variables in blank common. 
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Significant internal variables: 


FNT 

MXREG 

LIM1 

EB0T1 

EB0T2 

NGP3 

NQP1 

IFLG 


total number of source particles in the run, 
the number of geometry regions, 

the group number for the last primary 
particle group, 

lowest primary particle group energy, 
lowest secondary particle group energy, 
the group number for last secondary group, 

the group number for last primary group 
(negative if gamma-ray only problem), 

flag to indicate which fluence estimate is 
currently being output. 0-collision density, 
1-track length. 








Subroutine ENDRUN 


©■ 


(Collision density?) -Yes— 4~Set titles] 


1 Set titles .| 


(Calculate lowest neutron and gamma-ray energy, j 



(Calculate group widths, AE.J 



(Adjoint problem?)-Yes-~] Set AE = 1| 


IADJ >0 

No 


(Call CTV0L calculates 1/volume of each region. ( 


0 - 



































III. 


MODIFICATIONS FOR UNIVAC-1108 


In converting the M0RSE-SAMB0-M0RSEC code frorr the 
IBM 360 to the UNIVAC-1108, there were a considerable number of 
changes that had to be made. These changes were connected mostly 
with the local library functions of the IBM-360 and an attempt was made 
to not add local UNIVAC-1108 library functions. The only reduction in 
the code's capability is in the diagnostic module where the ability to 
look at a bit string and decide whether it was a floating point or integer 
number, convert it to hollerith and output it, does not exist in this version. 
However, the diagnostic module is mostly machine independent and still 
has the capability of writing out values from common and parts of blank 
common. Subroutines BNKHLP and HELPER are dummy routines. A 
diagnostic routine for the combinatorial geometry package has been added. 

Routines which are still machine dependent (involving coding 
other than word lengths, data statements and equivalences) are DATE 
and TIMER. Subroutine DATE calls ERTRAN and uses DEC0DE. Sub¬ 
routine TIMER c. Us CPUTIM that assumes that the time is returned in 
units of microseconds. 

Subroutine READSG is still machine dependent in its need to 
read each card twice and subroutine ERR0R has been added—its only 
function is to call EXIT. Subroutine INPUT now calls INPUT1 and 
INPUT2 with this change, permitting a better use of overlay. The 
first six variables in NXTRA from common AP0LL0 are used. 

The layout of blank common is given in Figs. 12-18 for the 
different modules of M0RSE. Definitions of variables in blank common 
are given in Table HI which defines the indexing scheme. Variables in 
the labelled commons are defined in Tables IV - XVI. 







A description of the energy indexing scheme in M0RSE- 
SAMB0-M0RSEC is given in Appendix A. 

3.1 Built-In Fluence Estimators 

There are two fluence estimators that are "built-in" the 
UNIVAC 1108 combinatorial geometry version of M0RSE. The addition 

of the collision density estimator to M0RSE has been discussed pre- 
( 1 ) 

viously v ' and the track length per unit volume estimator has been added 
to take advantage of a feature of the combinatorial geometry. The re¬ 
sulting fluence is averaged over geometry regions and if no regions are 
used only the spectrum averaged over the whole system is available. 

3. 2 Collision Density Fluence Estimator 

An estimate of the fluence averaged over a region specified 
by the geometry input may be obtained from the weight that has been 
summed over all collisions in the region by 



where 

cPy(EJ = fluence in region V at energy E., 

Ny. - number of particles of energy E i scatter' 

1 ing in region V, 

w = weight of the n th particle of energy E. 

n scattering in region V, 5 

NITS*NSTRT = total number of particles, 

VOly = volume region V, 


Z TV^V 


tota 1 cross section for groups E. in region 
V. 1 







3.3 


Track Length per Unit Volume Estimator 


Another estimate of the fluence averaged over a region may 
be obtained from the track length of the neutron or photon in that region. 
The track length calculated by the combinatorial geometry package is 
stored in a manner similar to the collision weights. The fluence esti¬ 
mate may be obtained by 


where 


N Vi 

„?x 


L w„ 


V E i> ■ NSWr*MTS*VOL„ 


= number of trajectories of energy E. in region V, 

w n = weight of particle on n th trajectory, and 

L = track length of n th trajectory in region V at energy 

n E. . 

1 


Subroutine ENDRUN is called to process the weights and 
track lengths to obtain estimates of the fluence. In an adjoint problem, 

the importance is calculated. The fluence estimate has units of particles/ 

2 2 
eV/cm /source and in an adjoint problem the units are importance/cm / 

source. The total fluence is integrated over energy from the primary 
and secondary energy groups to obtain an integral value for each region. 
The only input required for this routine is a data statement giving the 
medium number for each region; there is a limitation of only one medium 
^gion. In addition, one of the several options for inputing or calcu¬ 
lator-, .he volume of each region must be utilized. 

The major disadvantage with using these estimators is that 
no estimate of statistics is available. The use of this estimator is in¬ 
expensive, however (output routine is called only once), and is meant 
to be used as a supplemental estimator. 
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Starling 

Location 

Mnemonic: Variable 
Name' 

Length 


ENER 



VEL _ _ 

! 4+NMTG 


FS 

L0CWTS 

BFR 

) 


Current Weight 
Standards 

j 3+MGPREG 


Path 

| MGPRFG 


SPLIT and R.R. 
Counters 

J 8*MGPREG 

L0CFWL__ 

Initial Weight 
Standards 

J 3*MGPREG 


_Cun'ent.r\VL0_ 

initial FWLD 

J 2*MXREG 

L.0CEPR 

GW 1,0 

{ NMTG*MXREG 

L0CNSC . . 

EPRB 

} NMTG*MXREG 


Scattering 

Counters 

} 10*NMTG*MXREG 

L0CrSK . .- 

NSCA 

} MEDIA 


FISH 

} NMTC*MEDIA 


FSE 

| NMTG* MEDIA 

NGE0M. 

GMGN 

| NMTG+MEDIA 

NGLART_J 

Geometry Data 

j See Figs. 13-15 for 
’ details. 

NS1GL 

Permanent Cross 
Sections 

| See Fig. 16 for details. 

NMXSEC 

Temporary 

Cross 

Particle 

Bank 

J 14-rNMOST see Fig. 17 
| for details. 

NLAST 

Fission 

Bank 

\ 7*NMOST IF MFISTP>0 

J see Fig. 18 for details. 


User Aren 

} NLEFT see Fig. 19 for 
j details. 


Fig. 12. General layout of blank common. 
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Starting 

Location 


information 


NGE0M=NADD 

KMA 


KFI J D 

KLCR 


KNBD 


K30R 


KNSR 

KV0L 

NGIAST 


Length of geometry array 


MA ' 

| Integer array j 

FPD 

Floating point array 
L0CREG 

Indices to correlate MA array 
data wi'n code zone data ■ 

NUlto0D 

Number of bodies for each 
code zone 

IR0R 

Indices to correlate input zone 
to code zone i 

——t- 1 -\ 

MRIZ j 

Indices to correlate M0RSE 
region to input zone 

MRCZ 

Indices to correlate M0RSE 
region to code zone 

NMIZ 

Indices to correlate M0RSE 
media to input zone I 

NMCZ 

Indices to correlate M0RSE 
media to code zone 

KRlJ i 

Indices to correlate first code 
zone to ‘nput zones 

KR2 

Indices to correlate last code 
zone to input zone 

NS0R 

lndr.es of code zones m which 
source particles have been found 

VN0R 

Volume of each M0RSF region 


1 

vrjviA 

LFPD 

NUMR 


Fig. 13. Layout of combinatorial geometry data in bla-.k common. 











Figure 15, 







Position in 

Blank Common Information Stored Size Description 
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Data on last body in the zone 4 zones are searched. 







Position in 









r 


Location _ 

1RSG ®NLAST 

= NGLAST+1 

INGS 

INSG 

IDEL 

INNN 

IGGG 

ISTART 

ISCC0G 

INAB0G 

IGAB0G 

IFP0RG 

INUS 

INGN 

INGNP 

IFNGT> 

IFSP0G 

IDSG0G 

IPRBNG 

IPRBGG 

ISCANG 

ISCAGG 

ISP0RG + ISTART 
If ISTAT > 0 ISP0RT 


NLAST or NNIC 
NNIC 


_ Information 

List of Mixing Table 


Index to E gg (Primary) 

Index to Eg g (Secondary) 

List of Element I. D. Numbers 


Number of Downscatters for each 
primary group 


Number of Downscatters for each 
secondary group _ 


V*T ~ 

V £ t _ 

y V S T _ 

_ 

E^/E.j, not used at present 

E not used at present 


g -g 

Repeat for next medium 


l ; l Coefficient Primary 
Pj Coefficient Secondary 
Repeat for P L Coefficient 
R epeat for next medium 
Index to point cross sections for 
each multigroup boundary 


Size _ 

3*NMIX 

NGP 

NGG 

NELEM*NC0EF 
if IXTAPE •> 0 

NGP 

NGG 

NTG 

NTG 

NTG 

NGP 

NTG 

NTG 

NGG 

NGG* NGP 
NGP*NGG 
NDSNGP 
NDSNGG 
NDSNGP* NSCT 
NDSNGG* NSCT 
NDSNGP* NSCT 
NDSNGG* NSCT 
ISP0RG 
NDSNGP 
NDSNGG 

INFP0G*(NC0EF-2) 
INFP0G*(NC0EF-1) 
(IGQPT+1)*NMED if 
I06RT>O; otherwise 
zero 
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Location 

Information 

Size 

NXSECT(l) 

Point cross 
sections for 
Medium 1 

total 

scattering 

NPT (1)* (NEGPS-1)* 

NXSECT(J) 

Point cross 
sections for 
Medium J 

total 

scattering 

NXPM 

NPT(J)* (NEGPS-1)* 

NLAST 


" S f 

NXPM 


Fig. 16. Layout of permanent cross sections in blank common. 
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Location 


Variable 


■;N0 = NSIGL 


NN0 + 14 



Variables are from 
NUTR0N common, 
see Table VI for 
definitions. 


Repeat for Particle 2 


NN0 + 14*NM0ST = NLAST 

= NFISBN IF MFISTP > 0 


Fig. 17. Layout of particle bank in blank common. 












Location 

Variable 

NFISBN 

X 


Y 


Z 


WATEF 


AGE 


IG 


NAMEX 

NFISBN + 7 

Repeat for particle 2 i 

NFISBN + 7*NM0ST = NLAST 



Fig. 18. Layout of fission bank in blank common. 


L 
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Location Labels 

Mnemonic 
Variable Name 

Length 

NLAST* + 1 " — 



L0CRSP 

LABELS 

20*NRESP 


RESP 

NRESP*NMTG 

L0CXD 

EXTR 

NEX*NMTG 


XD 

YD 

ZD 

RAD 

T0 

FACT 

6*ND 

L0CIB 

EXTR 

NEXND*ND 


IB 



EP 

3*NE 

L0CC0 

DELE 


LOOT 

C0S 

NA 

L0CUD 

T 

DELT 

2*ND*NT 


UD 



SUD 

3*ND*NRESP 

L0CSD 

SUD2 



SD 



SSD 

SSD2 

3*ND*NRESP 

LOCOE 

SUD & SSD Units 

20 


NOTE: NLAST must not be changed after SC0RIN is called. 






















Indexing of Random Walk Blank Common Arrays 
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WRSU(IG,NREG) 


TABLE in. (Cont'd.) 
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_TABLE IV. (Cont'd.) 

Mnemonic Variable ' Location of Blank Common 

Name [BC(I) or NC(I)] Description 
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Response function. 


TABLE IV. (Cont’d.) 



75 



TABLE V. 


Definition of Variables in AP0LL0 Common 


Variable Definition 


AGSTRT 

DDF 

DEADWT(5) 

ETA 

ET/.TH 

ETAUSD 


Input starting age of source particle. 

Starting particle weight as determined in S0RIN. 

The summed weights of the particles at death. The 
four deaths are: Russian roulette, escape, energy, 
and age limit. DEADWT(5) is unused. 

Mean-free-path between collisions. 

Distance in cm to the next collision if the particle 
does not encounter a change in total cross section. 

Flight path in m. f. p. that has been used since the 
last event. 


UINP, VINP, WINP 
WTSTRT 

XSTRT, YSTRT, 
ZSTRT 

TCUT 

XTRA(IO) 

10, II 
MEDIA 
TADJM 
ISBIAS 

IS0UR 

ITERS 

ITIME 

ITSTR 


Input direction cosines for source particle. 

Input starting weight. 

Input starting coordinates for source particle. 

Age limit at which particles are retired. 

Not used. 

Output and input logical units. 

Number of media for which there are cross sections. 
Switch indicating an adjoint problem if > 0. 

Switch indicating that source energy distribution is 
to be biased if >0. 

Input source energy group if >0; otherwise, S0RIN 
is called to read input spectrum. 

Number of batches still to be processed in the run. 
Not. used. 

Switch indicating that secondary fissions are to be 
the source for the next batch if > 0. 
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TABLE V. (Cont'd.) 


Variable 


Definition 


L0CWTS 

L0CFWL 

L0CEPR 

L0CNSC 

L0CFSN 

MAXGP 




MEDALB 

MGPREG 

MXREG 

NALB 

NDEAD(5) 

NEWNM 

NGE0M 

NGPQT1* 

NGPQT2* 

NGPQT3* 


Starting location in blank common of the weight standards 
and other arrays MGPREG long. 

Starting location in blank common of the fission weights. 

Starting location in blank common of the energy-biasing 
parameters. 

Starting location in blank common of the scattering 
counters. 

Starting location in blank common of the fission and 
gamma-generation probabilities for ea-.h medium and 
group. 

Maximum number of energy groups for which there are 
weight standards or path-length stretching parameters. 

The elapsed clock time at which the problem is termin¬ 
ated. 

Medium number for the albedo medium. 

Product of number of weight standard groups (MAXGP) 
and regions (MXREG). 

Maximum number of regions in the system. 

An index indicating that an albedo scattering has occurred 
if >0. 

Number of deaths of each type (see DEADWT). 

Name of the last particle in the bank. 

Location of first cell of geometry data storage in blank 
connnc 

The lowest c-nergy group (largest group number) for 
which primary particles are to be followed. 

The number of primary particle groups. 

The lowest energy group (largest group number) for 
which any particle is to be followed. 


See Appendix A for diagram of energy group structure. 
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TABLE V. (Cont'd.) 


Variable 


Definition 


NGPQTG 

Number of energy groups of secondary particles to be 
followed. 

NGPQTN* 

Number of energy groups of primary particles to be 
followed. 

NITS 

Number of batches per run. 

NKCALC 

The first batch to be used for a k calculation. If 0, k 
is not calculated. 

NK3LL 

An index to indicate that Russian roulette is to be 
played if > 0. 

NLAST 

The last cell in blank common that was used by the 
cross-section storage or is set aside for banking. 

NMEM 

The location of the next particle in the bank to be pro - 
cessed. 

NMGP* 

The number of primary particle groups for which there 
are cross sections. 

NM0ST 

The maximum number of particles that the bank can hold. 

NMTG* 

The total number of energy groups (both primary and 
secondary'; for which there are cross sections. 

N0LEAK 

An index which indicates that nonleakage path-length 
selection is to be used if >0. 

N0RMF 

An index to indicate that the fission parameters are to 
be renormalized if >0. 

NPAST 

An index to indicate that the exponential transform is 
to be used if >0. 

NPSCL(13) 

An array of counters of events for each batch: 

1. sources generated 

2. splittings occurring 

3. fissions occurring 

4. gamma rays generated 


* 

See Appendix A for diagram of energy group structure. 
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TABLE V. (Cont’d.) 


Variable 


Definition 


NQUIT 

NSIGL 

NS0UR 

NSPLT 

NXTRA(l-5) 

NXTRA(6) 

NXTRA(7-10) 


5. real collisions 

6. albedo scatterings 

7. boundary crossings 

8. escapes 

9. energy cutoffs 

10. time cutoffs 

11. Russian roulette kills 

12. Russian roulette survivors 

13. gamma rays not generated because bank was full. 
Number of runs still to be processed. 

Starting location of the bank in blank common. 

An index input to indicate that fissions are to be the 
source for future batches. 

An index to indicate that splitting is to be considered 
if >0. 

Used to separate INPUT1 and INPUT2. 

NGPREG calculated in M0RSE and used in G0MST as 
an index for track length estimation. 

Not used. 




TABLE VI. 


Definition of Variables in NUTR0N Common 


Variable 

Definition 

NAME 

Particle's first name. 

NAMEX 

Particle's family name. (Note that particles do not 
marry.) 

IJ 

Current energy group index. 

IG0 

Previous energy group index. 

NMED 

Medium number at current location. 

MED0LD 

Medium number at previous location. 

NREG 

Region number at current location. 

U, V, w 

Current direction cosines. 

U0LD, V0LD, 
W0LD 

Previous direction cosines. 

X, Y, Z 

Current location. 

WATE 

Current weight. 

0LDWT 

Weight at previous collision. 

WTBC 

Weight just before curreni collision. 

BLZNT 

Current block and zone number (packed) for 05R geo¬ 
metry, and current zone number IR for combinatorial 
geometry. 

BLZ0N 

Previuus block and zone number (packed) for 05R geo¬ 
metry, and previous zone number IR for combinatorial 
geometry. 

AGE 

Current age. 

0LDAGE 

Previous age. 








TABLE VD. 


Definitions of Variables in USER Common 


AGSTRT 

WTSTRT 

XSTRT 

YSTRT 

ZSTRT 

DFF 

EB0TN 

EB0TG 

TCUT 

10 

II 

IADJM 

NGPQT1 

NGPQT 

NGPQT3 

NGPQTG 

NGPQTN 

NITS 

NLAST 

NLEFT 

NMGP 

NMTG 


Initial chronological age to be assigned to source 
particles. 

Initial weight to be assigned to source particles. 

Initial x position to be assigned to source particles. 
Initial y position to be assigned to source particles. 
Initial z position to be assigned to source particles. 
Normalization for adjoint problems. 

Lower energy boundary (eV) of last neutron group 
(group NMGP). 

Lower energy boundary (eV) of last gamma-ray group 
(group NMTG). 

Chronological age limit. 

Logical unit for output. 

Logical unit fcr input. 

Adjoint switch (>0 for adjoint problem). 

Problem-dependent energy group limits (see Appendix 
A). 

Group number of lowest energy gamma-ray group to be 
treated. 

Group number of lowest energy neutron group to be 
treated. 

Number of batches to be run. 

Last cell in blank common used by random walk package. 
Number of cells in blank common available to user. 
Number of primary particle energy groups. 

Total number of energy groups. 





TABLE VIII. 


Definitions of Variables in Common G0ML0C 


Variable 


Definition 


KMA 

Starting location for the array MA containing integer 
data for each code zone. 

KFPD 

Starting location for the array FPD containing real data 
for each code zone. 

KLCR 

Starting location for the array L0CREG(I) which contains 
the starting location in the MA array for the I** 1 code 
zone data. 

KNBD 

Starting location for the array NUMB0D(I) which con¬ 
tains the number of bodies for the I th code zone. 

KI0R 

Starting location for the array IR0R(I) which contains 
the index of the corresponding input zone for the I th 
code zone. 

KRIZ 

Starting location for the array MIRZ(I) which contains 
the index of the M0RSE region corresponding to the I th 
geometry input zone. 

KRCZ 

Starting location for the array MRCZ(I) which contains 
the index of the M0RSE region corresponding to the I th 
geometry code zone. 

KMIZ 

Starting location for the array MMIZ(I) which contains 
the index of the M0RSE media corresponding to the I^ 1 
geometry input zone. 

KMCZ 

Starting location for the array MMCZ(I) which contains 
the index of the M0RSE media corresponding to the I^ 1 
code zone. 

KKR1 

Starting location for the array KR1(L) contains the first 
code zone which was made from the L tJl input zone. 

KKR2 

Starting location for the array KR2(L) contains the last 
code zone which was made from the L th input zone. 

KNSR 

Starting location for array NS0R which contains the code 
zones in which source particles have been found. 


82 







TABLE Vm. (Cont’d.) 


Variable 

Definition 

KV0L 

Starting location for the array VN0R(I) which contains 
the volume for M0RSE region (I). 

NADD 

Starting location for the geometry data length and 
changed in J0MIN to the total number of words re¬ 
quired for geometry data. 

LDATA 

Length of the integer data in the MA array excluding 
the words set aside for zone search information. 

LTMA 

Total length of the MA array. 

LFPD 

Length of the FPD array 

NUMR 

Number of code produced zones. 

IRTRJ 

Number of input zones. 

NUMB 

Number of bodies. 

NIR 

Number of M0RSE geometry regions. 
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TABLE IX. 


Definitions of Variables in Common L0CSIG 


Variable 

Definition 

ISTART 

Starting location for the total cross-section vector for 
the first medium. 

ISCC0G 

Starting location . the scattering cross-section veccor 
for the first medium. 

INAB0G 

Starting location for the non-absorption vector for the 
first medium. 

IGAB0G 

Starting location for the gamma-ray production vector 
for the first medium. 

IFP0RG 

Starting location for vZf, the fission neutron production 
vector for the first medium. 

IFNGP 

Starting location for the primary-secondary transfer 
probability matrix. 

IFSP0G 

Starting location of the primary downscatter probability 
matrix. 

IDSG0G 

Starting location of the secondary downscatter prob¬ 
ability matrix. 

IPRBNG 

Starting location of the primary scattering angle prob¬ 
ability matrix. 

IPRBGG 

Starting location of the secondary scattering angle prob¬ 
ability matrix. 

ISCANG 

Starting location of the primary scattering angle matrix. 

ISCAGG 

Starting location of the secondary scattering angle 
matrix. 

ISP0RG 

Size of storage needed for each medium, not including 
Legendre coefficients. 



TABLE IX. (Cont'd.) 


Variable 


Definition 


ISP0RT* 

INPBUF 

ISIG0G 

INFP0G* 

IABS0G 

IT0TSG* 

NGP 

NDS 

NGG 

NDSG 

INGP 

INDS 

NMED 

NELEM 


Starting location for temporary storage of downscatter 
matrix. 

Starting location for temporary storage of the Pq table. 

Starting location for temporary storage of total cross 
section for element 1. 

Starting location for temporary storage of vLf for 
element 1. 

Starting location for temporary storage of downscatter 
matrix for P^ coefficients (primary groups, element 1). 

Total storage required for temporary storage. 

The number of primary groups to be treated. 

Number of downscatters for NGP (usually equal to NCP). 
Number of secondary groups to be treated. 

Number of downscatters for NGG (usually equal to NGG). 

Number of groups for which cross sections are to be 
input. 

Number of downscatters for the INGP groups. 

Number of media for which cross sections are to be 
stored--should be same as MEDIA as read on Card B 
of M0RSE input. 

Number of elements for which cross sections are to be 
read. 


If Legendre coefficients are to be restored, then: 

INFP0G - Redefined by JNPUT as number of locations required 
for each coefficient (both primary and secondary). 

IT0TSG - Redefined by JNPUT as total storage required for all 
coefficients for each medium. 

ISP0RT - Redefined by JNPUT as starting location of Pj co¬ 
efficient for primary groups for medium 1. 


85 










TABLE IX. (Cont'd.) 


Variable Definition 

NMIX Number of elements times density operations to be per¬ 

formed. 

NC0EF Number of coefficients, including Pq . 

NSCT Number of discrete angles (usually NC 0 EF / 2 i n t e g ra j)- 

NTS Number of downscatters for combined primary and 

secondary groups (usually equal to NTG). 

NTG Total number of groups (primary + secondary) = NGP 

f NGG. 

NDSNGP The number of locations needed for tne tiownscatter 

matrix for the primary particles. 

NDSNGG The number of locations needed for the downscatter 

matrix for the seconc \ry particles. 

IADJ Same as IADJM. 

NME Indicator for stripping gamma rays from a coupled 

neutron gamma-ray cross-sectior set--set equal to 
number of neutrons groups + 1. 

L0C Same as L0CEPR. 

INGS Starting location of the indices f< r starting location of 

the downscatter vector for each group for primar y 
particles. 

INSG Starting locat’- of the indices for starting location of 

tne downscattcx vector for each group for secondary 
particles. 

II, 10 Input and output logical unit numbers. 

KKK A running index of the number of cross sections that 

have already been read in (used in checking the element 
numbers obtained from tape). 

IXTAPE Logical tape number of the multigroup cross-section 

tape if > 0, or logical tape number of the processed 
cross-section tape if < 0. 

IDEL Starting location for element identifiers which determine 

the element cross sections to be read from tape. 
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Variable 


Definition 


ITEML 

ITEMG 

IRSG 

IRDSG 

ISTR 

IPRIN 

IFMU 

IM0M 

IDTF 

1STAT 

IPUN 

NUS 

NGN 

IHT 

INUS 

INUSN 

INGN 


Amount of storage for primary and secondary group 
downscatters per element. 

Starting location for temporary storage of downscatter 
matrix for Pl coefficients (secondary groups) for 
element 1. 

Starting location of the mixing parameters. 

Switch to print the cross sections and to test the card 
sequence as they are read if >0 (test card sequence 
only if = 0, and does neither if <0). 

Switch to print cross sections as they are stored if 

> 0 . 

Switch to print angles and probabilities if >0. 

Switch to print intermediate results of n's calculation 
if >0. 

Switch to print moments of angular distribution if > 0. 

Switch to signal that input format is DTF-IV format if 
>0; otherwise, ANISN format is assumed. 

Flag to restore Legendre coefficients for next-flight 
estimates if >0. 

Switch to print results of bad Legendre coefficients if 

> 0 . 


Number of groups of upscatter. 

Not use A presently. 

Location of L-j* in the cross-section table, currently 
set to 3 in XSEC. 

Starting location for the upscattering vector for the 
first medium. 

Number of impossible first moments found in the cross 
sections. 

Starting location for the photoneutron production vector 
for the first medium (not used at present—same as 
INUS). 
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TABLE IX. (Cont'd.) 

Defiiution 

Starting location for the secondary-primary transfer 
probability matrix (not used at present--same as INUS), 

Starting location for the array of the number of down- 
scatter groups for each primary group. 

Starting location for the array of the number of down- 
scatter groups for each secondary group. 


Variable 

INGNP 

INNN 

IGGG 








TABLE X. 


Definitions of Variables in Common GTSC1 


Variable 

Definition 

NPT(16) 

Number of points per supergroup for each medium. 

(Only 16 media are allowed.) Redefined at end of 

GTSCT as the total storage required for each cross 
section for each medium. 

NXSECT(17) 

Starting location of point cross sections for each 
medium. 

NDSGP 

Total storage required per supergroup for all media. 

11, 10 

Input and output logical unit numbers. 

I06RT 

Logical tape number of 06R point cross-section tape. 

IGQPT 

Last M0RSE multigroup for which point cross sections 
will be used (=£ NMGF'. 

NEGPS 

Number of supergroup boundaries. 

ESPD(IOO) 

Energy boundaries (in eV) of the supergroups (only 100 
supergroups are allowed). 

NNIC 

Starting location of the point cross sections. 

NXPM 

Number of cross sections per medium 

= 1 if total only 

= 2 if total + scattering 

= 3 if total, scattering, and y*fission. 

NINC 

Index for locating point cross sections. Used as flag in 
FISGEN to determine if transport process has begun. 


89 



TABLE XI. 


Definitions of Variables in Common PDET 


Variable 

Definition 

ND 

Number of detectors (2 1). 

NNE 

Number of neutron energy bins. 

NE 

Total number of energy bins. 

NT 

Number of time bins. 

NA 

Number of angle bins. 

NRESP 

Number of energy response functions 1). 

NEX 

Number of extra arrays of length NMTG to be set aside. 

NEXND 

Number of extra arrays of length ND to be set aside. 

NEND 

NE*ND. 

NDNR 

ND*NRESP. 

NTNR 

NT*NRESP. 

NTNE 

NT*NE. 

NANE 

NA*NE. 

NTNDNR 

NT*ND*NRESP. 

NTNEND 

NT*NE*ND. 

NANEND 

NA*NE*ND. 

L0CRSP 

Location of cell zero of response functions. 

L0CXD 

Location of cell zero of detector positions. 

L0CIB 

Location of cell zero of energy bins. 

L0CC0 

Location of cell zero of angle (cosine) bins. 

L0CT 

Location of cell zero of time bins. 

L0CUD 

Location of cell zero of array UD. 

L0CSD 

Location of cell zero of array SD. 

L0CQE 

Location of cell zero of array QE. 

L0CQT 

Location of cell zero of array QT. 
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TABLE XI. (Cont'd.) 


Variable Definition 

L0CQTE Location of cell zero of array QTE. 

L0CQAE Location of cell zero of array QAE. 

LMAX Last cell used in blank common. 

EFIRST Upper energy limit of first energy bin. 

EGT0P Upper energy limit of first gamma-ray bin. 








TABLE Xn. 


Definitions of Variables in Common FISBNK 


Variable Definition 


MFISTP Index to indicate that fissions are to be treated if > 0. 

NFISBM Location in blank common of cell zero of the fission 

bank. 

NFISH Number of fissions produced in the previous batch 

(generation). 

FT0TL Total fission weight from all collisions in this batch 

(generation). 

FWATE Total weight of fission neutrons stored in bank in this 

batch (generation). 

WATEF Weight of fission neutron stored in bank. 
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TABLE Xm. 


definitions of Variables in Common GE0MC 
as foirnd in 05R Geometry 


Variable Definition 


X2, Y2, Z2 Coordinates at tentative end-of-flight or if the trajectory 
is in an internal void; X2, Y2, Z2 are the direction 
cosines of the trajectory. 

XI, Yl, ZI Starting coordinates for the particle. 

ETA Number of mean free paths to be traversed if flight goes 

to X2, Y2, Z2. 

ETA USD Number of mean free paths actually traversed after the 

call to GE0M. 

IBLZ An index to the medium number for the special geometry 

packages. For GENERAL GE0M, IBLZ is a packed 
work giving the block and zone of the end of flight. 

IBZN A dummy variable. 

MARK A flag set by GE0M indicating the results of the trajec¬ 

tory calculation 

= 1 for completed flight 
= 0 for boundary crossing 
= -1 for escape 

= -2 for entering an internal void. 

NMED Medium number at end of flight or of medium about to 

be entered at a boundary crossing. 

NREG Region number at end of flight; not set at boundary 

crossings. 
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TABLE XIV. 


Definitions of Variables in Common PAREM 
as found in Combinatorial Geometry 


Variable 


Definition 


XB(3) Coordinates of the starting point of the present path. 

WB(3) Direction cosines of particle trajectory. Equal to U, 

V, and W. 

IR Combinatorial zone of present particle position. 

WP(3) Temporary storage of WB(3). 

XP(31 Temporary storage of XB(3). 

IDBG Set non-zero to initialize a debug printout. 

IRPRIM Next region to be entered after a call of Gl. 

NASC Body number of last calculated intersection. Set negative 

to indicate source or collision point not on a body sur¬ 
face. 

LSURF Surface of body NASC where intersection occurred. 

Positive if particle is entering the body and negative 
when exiting. 

NB0 Body number and a sign used to define zones. Input in 

zone description as positive when zone is contained in 
body and as negative if zone is outside body. 

LRI Entry surface calculated in GG. 

LR0 Exit surface calculated in GG. 

RIN Distance to entry calculated in GG. 

R0UT Distance to exit calculated in GG. 

KL00P Trajectory index of present path incremented in Gl. 

L00P Index of last trajectory calculated for body NB0. If 

L00P is equal to KL00P, GG returns immediately with 
old values of RIN, R0UT, LRI, and LR0. 

ITYPE Body type of body NB0 (indicates B0X, SPH, etc.). 
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TABLE XIV. (Cont’d.) 


Variable 

Definition 

PINF 

Machine infinity. 

N0A 

Not used. 

DIST 

Distance from XB(3) to present point. 







TABLE XV. 


Definitions of Variables in Common 0RGI 


Variable 

Definition 

M.'L.X 

Set 1 in G1 if trajectory end point is reached before 
next intersection. Otherwise set to 0. 

DISTO 

Distance from point XB(3) to next scattering point. 

Used in G1 to avoid calculating the next zone if a 
scattering event occurs before the intersection. 

NMEDG 

Zone number IR from a L00KZ call. Stored in BLZNT 
by MS0UR. 


Note: Variabl names are not the same in all routines. 








TABLE XVI. 


Definitions of Random Walk Variables in Blank Common 


Mnemonic 

Variable Definition 

Name 


ENER(IG) Upper energy boundary of group IG (in eV). 

VEL(IG) Velocity corresponding to the mean energy for neutron 

groups and the speed of light for gamma-ray groups 
(in cm/sec). 

FS(IG) Unbiased source spectrum—unnormalized fraction of 

source particles in each energy group--transformed to 
c. d. f. by S0RIN. 

BFS(IG) Biased source spectrum--relative importance of each 

energy group—transformed to biased c. d. f. by S0RIN. 

WTHI(IG, Weight above which splitting is performed (vs. group 

NREG) and region). 

WTL0(IG, Weight below which Russian roulette is performed (vs. 
NREG) group and region). 

WTAV(IG, Weight tc be assigned Russian roulette survivors (vs. 
NREG) group and region). 

PATH(IG, 

NREG) Exponential transform parameters (vs. group and region). 
NSPL(IG, 

NREG) Splitting counter (vs. group and region). 

WSPL(IG, 

NREG) Weight equivalent to NSPL. 

N0SP(IG, Counter for full bank when splitting was requested (vs. 

NREG) group and region). 

WN0S(IG, 

NREG) Weight equivalent to N0SP. 

RRKL(IG, 

NREG) Russian roulette death counter (vs. group and region). 








TABLE XVI. (Cont'd.) 


Mnemonic 

Variable Definition 

Name 


WRKL(IG, 

NREG) 

RRSU(IG, 

NREG) 

WRSU(IG, 

NREG) 

INIWHI (IG, 
NREG) 

INIWL0(IG, 

NREG) 

INIWAV(IG, 

NREG) 

FWL0(NREG) 


INIFL0(NREG) 

GWL0(IG, 

NREG) 

EPRB(IG, 

NREG) 

NSCT(IG, 

NREG) 

WSCTdG, 

NREG) 

NALBdG, 

NREG) 

WALBdG, 

NREG) 

NFIZdG, 

NREG) 


Weight equivalent to RRKL. 

Russian roulette survival counter (vs. group and 
region). 

Weight equivalent to RRSU. 

Initial values of WTHI array. 

Initial values of WTL0 array. 

Initial values of WTAV array. 

Weights to be assigned to fission daughters (vs. 
region). 

Initial values of FWL0. 

Weights to be assigned to secondary particles (vs. 
group and region). 

Relative importance of energy groups after scatter¬ 
ing (vs. group and region). 

Number of real scatterings (vs. group and region). 

Weight equivalent to NSCT. 

Number of albedo scatterings (vs. group and 
region). 

Weight equivalent to NALB. 

Number of fissions (vs. group and region). 






TABLE iVI. (Cont'd.) 


Mnemonic 

Variable 

Name 

Definition 

WFIZ(IG, 

NREG) 

NGAM(IG, 

NREG) 

Weight equivalent to NFIZ. 

Number of secondary productions (vs. group and region). 

WGAM(IG, 
,NREG) k 

Weight equivalent to NGAM. 1 

NSCA(IMED) 

Scattering counter (vs. cross-section medium). 

FISH(IG, 

IMED) 

Probability of generating fission neutron (vs. group 
and medium). 

FSE(IG, 

IMkD) 

Source spectrum for fis$ion-induced neutrons for each 
group—input as frequency of group IG. 

GMGN(IG, 

IMED) 

Probability of generating secondary particle (vs. group 
and medium). 

j_ 



99 







IV. 


THE PICTURE PROGRAM 


This program is used to obtain printer plots of two-dimensional 
slices through the three-dimensional geometry described by the combi¬ 
natorial geometry (CG) package. The 05R geometry version* 1 ^ has 
been modified for use with CG and some minor changes to the program 
have been made. It is now possible to obtain a picture the full width of 
a computer page and the logic for determining the medium and/or region 
number has been changed so that a test of the geometry ray tracing cap¬ 
ability is also made. That is. instead of just testing the routine L00KZ, 
this modification forces the geometry package to execute the routines 
and logic used during the random walk process. This program has been 
very useful in debugging geometry input. It is suggested tjiat any com¬ 
plex geometry be tried in this program before using it in the M0RSE 
code. A sample problem is provided. Besides the geometry package, 
a main program, PICTURE, and two routines, MESH and PRINT, are re¬ 
quired. A discussion of these routines is included and a definition of 
the variables in common, PICT, is given in Table XVII. 


Preceding page blank 
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4.1 


Program PICTURE 

PICTURE is the executive routine for the PICTURE program 
and reads in the input, calculates the coordinates of the picture to be 
plotted and controls the calls to other routines. There are several dif¬ 
ferent ways in which a two-dimensional slice through the geometry may 
be obtained. These different options are discussed in Appendix B with 
discussions of the input. The characters to be printed for corresponding 
media or regions may be changed by altering the values in ATABLE as 
given in the data statement. 

Subroutines called: J0MIN, PRINT 
Commons required: PICT 

Variables required: Several input cards are read. 

Variables changed: All variables in common PICT. 

Significant internal variables: 

INT - input logical unit, 

I0T - output logical unit, 

NADD(l) - first location in blank common for storage 
of geometry data, 

flag to indicate which set of input options 
was used to define the two-dimensional 
slice. 


NCK 







Program PICTURE 


d> 

€> 


I Set input/output devices to 5 and 6, j 


(l/&. Sym bol option. | 
Read symbol tab l e if required") 


---i___ 

[ Call J0MIN to input geometry data.) 


V_ 


J 

_Z=J 

11/0: Options, number of pictures for this 

1 geometry, media or region title. 


fl/P: picture coordinates upper left and 
lower right corner. _ 


ft: direction numbers for U and v axes! 

1 - J 

fr/0: interval number and size.) 



fc print [ 


Q)*~ WCT>0 fls there another gcometryTnput ?) 
©- 1 
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4.2 


Su broutine PRINT (KXX, KYY, ATABLE) 


This routine controls the printing of the picture. First it is 
determined if the picture is to be more than one ’’page” (130 columns) 
wide and then for each line of the picture on the first page the infor¬ 
mation to be printed is determined by calling MESH and printed with a 
130A1 format. The next lines are then calculated and printed until that 
"page” is finished. Note that a page refers to width not length. Thus, 
as much detail may be obtained as necessary in both directions by piec¬ 
ing together the output. 

Called from: PICTURE 

Subroutine called: MESH 
Commons required: PICT 
Variables required: 


KXX 

number of intervals in the U direction 
(direction of paper movement through 
the printer), 

KYY 

number of intervals in the V direction 
(line), 

ATABLE 

table of characters to be printed. 


Significant interval variable: 

NPAGES - number of subpictures required to cover 
the width of the total picture, 

I0T - output logical unit, 

NV - number of characters per line (characters/ 

page width). 
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Subroutine PRINT 
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4. 3 Subroutine MESH (XS, YS, ZS, NV) 

Subroutine MESH is used by the PICTURE package to set up 
one line of print in the array NST0R. Both for efficiency and to debug 
the combinatorial geometry package, this version has been modified to 
work exactly like particle tracking. L00KZ is first called to deter¬ 
mine the zone of the first grid point. A trajectory to the last grid point 
is then initialized, and successive calls to G1 "track” a particle to the 
last point, setting the region of each grid point in NST0R. By setting 
IRG negative, .°ro, or positive, either NREG, IR, or NMED will be 
stored in the print array NST0R. 

Called from: PRINT 

Subroutines called: L00KZ, G1 

Commons required: PICT, G0ML0C, PAREM, 0RGI, blank. 

Variables required: 

XS, YS, ZS - coordinates of first grid point, 

NV - number of grid points, 

DELV - distance between grid points, 

IRG - flag to print NREG, IR, or NMED if IRG 

is negative, zero, or positive. 

Variables changed: NST0R - print array. 

Significant internal variables: 

J - grid point index, 

IST0R - value to be stored in NST0R between 

successive boundary crossings. 
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Subroutine MESH 
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TABLE XVH. 


Definition of Variables in Common PICT 


Variable 

Definition 

DELU 

The increment in geometry units between lines in the 
picture in the U direction. 

DELV 

The increment in geometry units between lines in the 
picture in the V direction. 

XO, YO, ZO 

The coordinates in geometry units defining the first 
point (upper left hand corner) of the picture. 

XU, YU, ZU 

The length (in direction U) of the picture in geometry 
units. 

XV, YV, ZV 

The width (in direction V) of the picture in geometry 
units. 

NST0R(13O) 

An array used to store the medium, region, or zone 
number for one line of the picture. 

IRG 

A flag indicating that region, zone, or medium geo¬ 
metry parameter should be printed if IRG is negative, 
zero, or positive, respectively. 




V. M0RSE SAMPLE PROBLEMS 

5.1 Sample Problem Number 1 

(7) 

This sample problem is similar to that used previously 
which calculates neutrons and secondary gamma rays in infinite homo¬ 
geneous air. The source is a 12. 2- to 15-MeV neutron source. Two 
types of estimators are used in the analysis of the particle histories. 
First, a boundary-crossing estimator is used for five spherical shells. 
(This part of the problem is similar to the M0RSE sample problem.) 
Second, a collision density estimator is used in RELC0L for nine de¬ 
tector regions (detectors 6-14). By choosing the spherical shells as 
the midpoints of the regions, an internal check is provided. Also, by 
choosing the geometry regions to correspond to the collision density 
detectors, the results from SAMB0 can be compared with results from 
ENDRUN. With the inclusion of the combinatorial geometry package, 
an additional estimator based on track length/unit volume has been 
added with results also output in ENDRUN. Both the collision density 
and track length per unit volume estimators are obtained for the geo¬ 
metry regions without s. ecific coding. 

Since the collision density estimator in RELC0L scores 

2 

4ffR. <p there will be some differences between these estimates and 
those from the collision ..ansity estimator output in ENDRUN which 
2 

scores 4ir <R> <p. (The symbols ( > denote an average over a spatial 

region.) 

Table XVIII illustrates the comparison for one of several 
sets of answers as calculated by the various estimators. Descriptions 
of the user routines INSC0R, RELC0L, and BDRYX required for the 
sample problem are given. 







TABLE XVm. 


Comparison of Results of Sample Problem 1 


Detector No. 

Type and Location 
of Detector 

Result 

Response #1 Group 1 * 

4ttr 2 fluence fluence/eV 
top group 

3 

Boundary crossing 
at 450 meters. 

8.19(-2) 

2. 92(-8)*A 
= 8.19(-2) 

9 

Collision density for 
volume between 400- 
500 meters. 

7. 59(-2) 

2. 71(-8)*A 
= 7. 59(-2) 

Region 4 

Volume between 
400-500 meters 
(ENDRUN results). 

Real 

Collision 
Region 4 

Top Group 

1.11(-18)*A*B 
= 7. 76(-2) 


Volume between 
400-500 meters 
(ENDRUN results). 

Track 

Length 

Region 4 

Top Group 

1.18(-18)*A*B 
= 8. 36(-2) 


A = width of the top group 2. 8(+6) eV. 
B = Alt ^ for 4** 1 region = 2. 5(+10). 





5.1.1 Subroutine INSC0R - Sample Problem Number 1 


This user routine is used to calculate the volume of each of 
the detector regions (detectors 6-14) in the sample problem. The re¬ 
ciprocal of the volume is stored in the array in blank common set aside 
for FACT, a detector -dependent normalization. Because two types of 
estimators are used in the problem, normalization values of 1 are 
stored for the boundary detectors. The detector number NDC for the 
last boundary crossing detector (in this case, 5) must be input. 

Called from: SC0RIN 

Subroutines called: None 

Commons required: PDET, USER, DETCUT, blank. 

Variables required: 

NDC - the number of the last boundary crossing 

detector. Detector radii from blank 
common. 

Variables changed: 

Normalization array in blank common. 


Ill 







Subroutine INSC0R 




5.1.2 Subroutine RELC0L - Sample Problem Number 1 


This subroutine estimates the iluence from the density of 
2 

collisions. For comparison purposes, 4 ttR times the fluence is stored 
where R. is the radius to each collision site. Thus, the contribution is 

4»r 2 ^5£ , 

1 S T 

where 

R. -- X 2 + Y 2 + Z 2 , 
l ill* 

WTBC = particle weight before collision, 

Z = cross section for the group (IG0) of the incoming 
particle. 

A detector number from NDC + 1 to ND is chosen depending 
on the value of R i . 

Called from: BANKR(5) 

Subroutines called: FLUXST, NSIGTA 

Commons required: PDET, NUTR0N, USER, DETCUT 

Variables required: 

X, Y, Z, IG0, NMED, AGE, WTBC from NUTR0N common 
TS1G and detector radii. 

Significant internal variables:- 

R - radial distance to collision site, 

C0N - fluence estimate. 
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Subroutine RELC0L 
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5.1.3 


Subroutine BDRYX - Sample Problem Number 1 ! 

This subroutine is called whenever the partiple in the random 
walk encounters a change in geometry media. If the source-to-collision 
distance corresponds to a detector position, the reciprocal of the cosine 
of the angle fij'om the radius vector is used as a fluence estimate. 

Called from: BANKR(7) 

Subroutines called: ERR0R 

Coriimons required: USER, NUTR0N, DET 

Variables required: 

X, Y, Z, U, V, W, WATE (from common NUTR0N). 
Significant internal variables: 


R21 

Radial distance to boundary crossing, 

R2 

- 99% of R21, 

R22 

- 101% of R21, 

C0S 

cosine of angle between particle direction 
and radius vector, 

ABC0S 

absolute value of C0S, 

C0N 

fluence estimate. 


\ 


\ 
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Subroutine BDRYX 
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5.2 


Sample Problem Number 2 


This sample problemis used to illustrate the next event 

estimator. The neutron and secondary gamma ray spectra at several 
4 

point detectors located 10 cm from a polyethylene slab is calculated. 
The source is a beam of neutrons with a fission spectrum energy dis¬ 
tribution. The use of the point detector routines to estimate the energy 
dependent fluence is illustrated and in particular, for secondary gamma 
rays the contribution of the uncollided or first flight gamma rays are 
scored in subroutine SGAM. The estimates from real collisions are 

made in subroutine RELC0L. Results are recorded in Table XIX in 
2 

units of particles cm source neutron. 


TABLE XIX. 


Results of Last Flight Estimator 
Sample Problem 


Detector 

2 

Neutrons/cm 

Per 

Source Neutron 

2 

Gammas/cm 

Per 

Source Neutron 

1 

6. 24(-ll) 

2. 36(-ll) 

2 

6. 44(-10) 

4. 82(-ll) 

3 

1. 32(-9) 

4. 93(-ll) 

4 

1. 56(-9) 

4. 94(-ll) 

5 

6. 92(-l) 

2.15(-3) 
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5.3 


Sample Problem Number 3 


This sample problem serves two purposes by (1) demonstrat¬ 
ing the use of the PICTURE code, and (2) verifying the tracking capability 
of the combinatorial geometry package in a moderately complex three- 
dimensional system. In the PICTURE program, subroutine G1 is called 
from MESH in the exact manner that it is called from G0MST during a 
M0RSE calculation. As such, PICTURE works by tracking rays and 
thus can be used to debug a combinatorial geometry input. 

The sample problem itself is a picture of a tank. The tank 
model itself was constructed purely as an illustration of the combina¬ 
torial geometry and is in no way accurate or detailed. Notice that the 
use of both the 0R operator and the ARB body is demonstrated. 
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TABLE XX. (Cont’d.) 

























































































APPENDIX A 


ENERGY INDEXING SCHEME 

For most problems the user does not need to worry about the 
energy-indexing scheme used in the M0RSE Monte Carlo code. Energies 
corresponding to the group boundaries in the forward group structure 
are input and the code takes care of the rest. The analysis package 
SAMB0 performs the proper bookkeeping for labeling. A diagram for 
the energy group structure is given in Fig. Al. However, some clarifi¬ 
cation of the indexing scheme for several types of problems is advanta¬ 
geous. This clarification can perhaps best be made through the use of 
examples. Consider a five-group coupled neutron-gamma-ray problem 
with the following group boundaries: 

Group Energy (eV) 

1 15(+6) 

2 5 (+6) 

3 l(+6) 

EB0TN 3 (+3) 

4 10 (+6) 

5 2 (+6) 

EB0TG 4 (+5) 

There are several options which might be considered; namely, 
(1) a forward coupled problem, (2) a forward neutron-only problem, (3) 
a forward gamma-ray-only problem, (4) an adjoint coupled problem, (5) 
an adjoint neutron-only problem, and (6) an adjoint gamma-ray-only 
problem. Table AI gives the energies as they are indexed during a 
M0RSE run. Table All gives the values of variables NQTl, NQT2, and 
NQT3 as they appear in USER common for the same six problems. 

Preceding page blank 


with NMGP = 3, NMTG = 5, 
NGPQTN = 3, NGPQTG = 2 
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Forward 


Adjoint 


NGPQT1 


NGPQT2 


NGPQT3* 


1 


NMTG 

“ 


* 

NGPQTN 



- 


NMTG - NGPQTN 

NMGP 


- 




- 


NMTG - NMGP 

- 


I 

NMGP + NGPQTG 



i 


NMTG - NMGP - NGPQTG 

- 

' 

- 

NMTG 


1 


Fig. Al. Diagram of energy group structure. 


NGPQT3 


NGPQT2 


- NGPQT1 


-NGPQTO 


NGPQT3 = NGPQT1 for neutron only or gamma ray only problem. 

















I 






TABLE AI. (Cont'd.) 















APPENDIX B 


INPUT INSTRUCTIONS 


B. 1 M0RSE Input Instructions 

The input read 1 by subroutine INPUT is as follows: 

CARD A (20A4) 

Title card. (Any character other than a blank or alphameric in 
column one will terminate the job.) 


CARD B (1315) 
NSTRT 
NM0ST 


NITS 

NQUIT 

NGPQTN* 

NGPQTG* 

NMGP* 


number of particles per batch, 

maximum number of particles allowed for in the \ 
bank(s); may equal NSTRT if there is no splitting, 
fission, and secondary generation during ex¬ 
ecution. If bank size is exceeded by more than 
50 due to fission or secondary gamma ray gen¬ 
eration the job is terminated, 

number of batches, 

number of sets of NITS batches to be run without 
calling subroutine INPUT, 

number of neutron groups being analyzed, 
number of gamma-ray groups being analyzed, 

number of primary particle groups for which 
cross sections are stored; should be same as 
NGP (or the same as NGG when NGP = 0) on 
Card XB read by subroutine XSEC, 


Input and output logical unit members are defined in MAIN and ii Gl. 
Unit 16 is a temporary storage unit used by J0MIN. 

* 

See Table BI for sample input. 
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NMTG* 

NC0LTP 

IADJM 

MAXTiM 

MEDIA 

MEDALB 


total number of groups for which cross sections 
are stored; should be same $s NGP + NGG as 
read on Card XB read by subroutine XSEC, 

set greater than zero if a collision tape is de¬ 
sired; the collision tape- is written by the user 
routine BANKR, 

set greater than zero for an adjoint problem, 

maximum clock time in minutes allowed for the 
I problem to be on the computer, 

number of cross-section media; should agree | 
w'th NMED on Card XB read by subroutine XSEC, 

albedo scattering medium is absolute value oi 
MEDALB; if 

= 0, no albedo information to be read in, 

<0, albedo only problem—no cross sections 
are to be read, 

> 0, coupled albedo and transport problem. 


CARD C (415, 5E10. 5) 
IS0UR 


ISBIAS 

N0TUSD 

WTSTRT 

EB0TN 

EB0TG 


source energy group if > 0, if IS0UR £ 0, 

S0RIN is called for input of Cards El and E2, 

number of groups for which the source spectrum 
is to be defined. If IS0UR s 0, NGPFS s 2, 

\ no source energy biasing if set equal to zero; 
otherwise the source energy is to be biased, 
and Cards E2 are required, 

an unused yariable, 

weight assigned to each source particle, 

lower energy limit of lowest neutron group (eV) 
(group NMGP), 

lower energy limit of lowest gamma-ray group 
(eV) (group NM7G), : 


See Table BI for sample input. 






TCUT - age in sec at which particles are retired; if 
TCUT = 0, no time kill is performed, 

VELTH - velocity of group NMGP when NGPQTN > 0; i. e., 
thermal-neutron velocity (cm/sec). 

i 

CARD D (7E10.4) 

XSTRT ) \ 

YSTRT > - coordinates for source particles, 

^ ZSTRT ) 1 

1 AGSTRT ' - starting age for source particles, 

source particle direction cosines; if all are 
zero, isotropic directions are chosen. 

Source data on Cards C and D may be overridden by changes 
in subroutine S0URCE. 

CARDS El .'7E10.4) (omit if IS0UR on Card C > 0) 

NGPFS values of FS, where FS equal the unnormalized fraction 
of source particles in each group. 

CARDS E2 (7E10. 4) (omit if IS0UR >0 or if IS0UR ^ 0 and ISBIAS = 0) 

If ISBIAS >0, NGPFS values of BFS, the relative importance of 
a source in group I, are required. 

CARDS F (7E10. 4) 

NMTG values of ENER, the energies (in eV) at the upper edge of 
the energy group boundaries. (Note: The lower energies of 
groups NMGP and NMTG were read on Card C.) 

CARD G (215, 5X, 3611, 1311) (omit if NC0LTP on Card B £ 0) 

NHISTR - logical tape number for the first collision tape, 

NHISMX - the highest logical number that a collision tape 

j may be assigned, 


UINP ) 
VINP J - 
WINP ) 
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NBIND(J), J-l, 36 - an index to indicate the collision para¬ 

meters to be written on tape, 

NC0LLS(J), J=l, 13 - an index to indicate the types of collisions 
to be put on tape. (See Table BII for in¬ 
formation concerning NBIND’s and 
NC0LLS.) 


CARD H (Z12) 

RAND0M - starting random number. 


CARD I (715) 
NSPLT 
NKILL 

NPAST 

N0LEAK 

IEBIAS 

MXREG* 

MAXGP 


index indicating that splitting is allowed if > 0, 

index indicating that Russian roulette is allowed 
if >0, 

index indicating that exponential transform is in¬ 
voked if > 0 (subroutine DIREC required), 

index indicating that non-leakage is invoked if 

> 0 , 

index indicating that energy biasing is allowed 
if >0, 

number of regions described by geometry input 
('- r ill be set to one if ^ 0), 

group number of last group for which Russian 
roulette or splitting or exponential transform is 
to be performed. 


I: ENDRUN is used, a data array relating media numbers to region 
numbers must be given in a data statement in ENDRUN. 
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CARD J (615, 4E10. 5) (omit if NSPLT + NKILL + NPAST = 0) 

from energy group NGPl to energy group NGP2, 
inclusive, in steps of NDG and from region NRG1 
to NRG2, inclusive, in steps of NDRG, the follow¬ 
ing weight standards and path-stretching para¬ 
meters are assigned. If NGPl = 0, groups 1 to 
MAXGP will be used; if NRG1 = 0, regions 1 to 
MXREG will be used (both in steps of one). Us¬ 
ually NDG = 1 and NDRG = 1. 

WTHIH - weight above which splitting will occur, 

WTL0W - weight below which Russian roulette is played, 

WTAVE - weight given those particles surviving Russian 

roulette, 

PATH - path-length stretching par. meters for use in ex 

ponential transform (usually 0 £ PATH < 1). 

The above information is repeated until data for all groups and 
regions are input. 

End Cards J with negative value of NGPl (ex., -1 in columns 4 
and 5). 

CARDS K (7E10. 4) (omit if IEBIAS on Card I * 0) 

((EPR0B(IG, NREG), IG = 1, NMTG), NREG = 1, MXREG) 

Values of the relative energy importance of particles leaving a 
collision in region NREG. Input for each region must start on a 
new card. 


CARD L (415) 

NS0UR - set £ 0 for a fixed source problem; otherwise the 
source is from fissions generated in a previous 
batch, 

MFISTP - index for fission problem, if £ 0 no fissions are 
allowed, 

NKCALC - the number of the first batch to be included in the 
estimate of k; if * 0 no estimate of k is made, 


NGPl 

NDG 

NGP2 

NRG1 

NDRG 

NRG2 
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N0RMF 


the weight standards and fission weights are un¬ 
changed if * 0; otherwise fission weights will be 
multiplied, at the end of each batch, by the latest 
estimate of k and the weight otandards are mul¬ 
tiplied by the ratio of fission weights produced in 
previous batch to the average starting weight for 
the previous batch. For time-dependent decay¬ 
ing systems, N0RMF should be > 0. 


CARDS M (7E10.4) (omit if MFISTP on Card L * 0) 

(FWL0(I), 1=1, MXREG) values of the weight to be assigned to 
fission neutrons. 


CARDS N (7E10.4) (omit if MFISTP on Card L^O) 

(FSE(IG, IMED), IG = 1, NMGP), IMED = 1, MEDIA) the fraction 
of fission-induced source particles in group IG and medium IMED. 

Note: Input for each medium must start on a new card. 


CARDS 0 (7E10. 5) (omit if NGPQTN = 0 or NGPQTG = 0, i. e., include 
if coupled neutron-gamma-ray problem) 

((GWL0(IG, NREG), IG = 1, NMGP or NMTG - NMGP), NREG = 1, 
_ average number of secondary particles to be 
generated. NMGP groups are read for each 
region in a forward problem and NMVG - 
NMGP for an adjoint. Input for each region 
must start on a new card. 
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TABLE BI. 


Sample Group Input Numbers for Some 
Representative Problems* 


Case A - Neutron Only Cross Sections (22 groups) 

Case B - Gamma-Ray Onlv Cross Sections (18 groups) 

Case C - Neutron-Gamma-Ray-Coupled Cross Sections (22-18 groups) 


Input Variable 

Problem Type 

Case A 

Top 14 groups 

Case B 

Top 17 groups 

Case C 

Neutrons Only 

Top 14 groups 

Case C 

Gamma Rays Only 
Top 17 groups 

Case C 

Neutron-Gamma 
Top 14 Neutron, 
Top 17 Gamma 
groups 


M0RSE Input: 







NGPQTN 

14 

0 

14 

0 

14 


NGPQTG 

0 

17 

0 

17 

17 

CARD B, 

NMGP 

22 

18 

22 

18 

22 

Variables 

NMTG 

22 

13 

22 

18 

40 


NGP 

22 

18 

22 

0 

22 


NGG 

0 

0 

0 

18 

18 

CARD XB 

INGP 

22 

18 

40 

40 

40 

Variables 


For cross sections with full downscatter, NDS = NGP, NDSG = NGG, 
INDS = iNGP, and ITBL = number of downscatters + number of up- 
scatters + 3. Usually, ISGG = number of upscatters + 4. 
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TABLE BE. 


BANKR Arguments (NC0LLS) 


BANKR 

Argument 

Called 

from 

Location of Call in Walk 

-1 

M0RSE 

After call to INPUT--to set parameters 
for new problem 

-2 

M0RSE 

At the beginning of each batch of NSTRT 
particles 

-3 

M0RSE 

At the end of each batch of NSTRT 
particles 

-4 

M0RSE 

At the end of each set of NITS batches— 
a new problem is about to begin 

1 

MS0UR 

After a source event 

2 

TESTW 

After a splitting has occurred—com¬ 
mented in column 1 

3 

FPR0B 

After a fission has occurred 

4 

GST0RE 

After a secondary particle has been 
generated 

5 

M0RSE 

After a real collision has occurred— 
post-collision parameters are avail¬ 
able 

6 

M0RSE 

After an albedo collision has occurred— 
post-collision parameters are available 

7 

NXTC0L 

After a boundary crossing occurs (the 
track has encountered a new geometry 
medium other than the albedo or void 
media) 

8 

NXTC0L 

After an escape occurs (the geometry 
has encountered medium zero) 

9 

M0RSE 

After the post-collision energy group 
exceeds the maximum desired—com¬ 
mented in column 1 
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TABLE BD. (Cont'd.) 


BANKR 

Argument 

Called 

from 

Location of Call in Walk 

10 

M0RSE 

After the maximum chronological age 
has been exceeded—commented in 
column 1 

11 

TESTW 

After a Russian roulette kill occurs— 
commented in column 1 

12 

TESTW 

After a Russian roulette survival 
occurs—commented in column 1 

13 

GST0RE 

After a secondary particle has been 
generated but no room in the bank is 
available—commented in column 1. 
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1 


B. 2 05R Geometry Input Instructions 

Read by subroutine J0MIN and the specific input depends on the 
geometry packages used. 

B. 2.1 Spherical GE0M 
CARD GA (15, DIO. 5) 

MED - medium number interior to R(> 0), 

R - outer radius of sphere or spherical shell contain¬ 

ing MED. 

Repeat Card GA for all radii (s 20) in increasing order. End 
Card GA input with blank card. 


CARD GB (DIO. 5) 

R - region radius of sphere or spherical shell con¬ 

taining regions. Region numbers are assigned 
in consecutive order starting with 1, and R must 
be in increasing order. 

Repeat Card GB for all radii (* 20) in increasing order. End Card 
GB input with blank card. If no regions are desired, a blank card 
must be used to signal no region geometry. 

B. 2.2 Slab GE0M 

CARD GA (15, DIO. 5) 

MED - medium with Z as lower bound (>0), 

Z - lower limit of medium MED. 

Repeat Card GA for all boundaries with the last card containing 
MED = 0 and the boundary of the system. 


J- 


138 







CARD GB (DIO. 5) 

Z - lower limit of region boundary. Region numbers 

are assigned in consecutive order starting with 
1, and Z must be in increasing order. 

Repeat Card GB for all region boundaries. 

End Card GB input with a blank card. If no region geometry is 
desired, a blank card is required. 


CARD GC (4D10. 5) 

XL 

XU 

YL 

YU 


lower boundary of 
upper boundary of 
lower boundary of 
upper boundary of 


system in X direction, 
system in X direction, 
system in Y direction, 
system in Y direction. 


B. 2.3 Cylindrical GE0M 
CARD GA (15, 5X, A8) 

NREGIN - flag to indicate material media (=2) or both 
region and material media (=1), 

SEX - sex of programmer. 

CARD GB (E10. 5) 

R - radii of the cylindrical shells describing the 

material media in ascending order. 

Repeat Card GB until all radii have been input. 

End Card GB input with a blank card. 
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CARD GC (EiO. 5, 1215/815) 


H - upper height of medium M(I) (> 0), cyUnders as¬ 

sumed to start at H = 0, 

M(I) - media for the cylindrical shells for this height. 

Repeat Card GC until all height intervals have been input. 

End Card GC input with a blank card or if there are more than 12 
radial intervals, 2 blank cards. 


CARD GC (E10. 5) (omit if NREGIN = 2) 

RG - radii of the cylindrical shells describing the 

region geometry in ascending order. 

Repeat Card GD until all region geometry has been input. 
End Card GD with a blank card. 


CARD GE (E10. 5, 1215/815) (omit if NREGIN = 2) 

HG - upper height of region MG(I), 

MG (I) - region numbers for the cylindrical shells for this 

height. 

Repeat Card GE until all height intervals have been input. 

End Card GE input with a blank card or if there are more than 12 
radial intervals, 2 blank cards. 


B. 2.4 General GE0M 


CARD GA (15, 5X, A6, IX, A7) hollerith left adjusted 


NSTAT 

SEX 

STATUS 


flag to indicate material media only if 2 and both 
region and material media if 1, 

sex of the programmer (select one from MALE, 
FEMALE, or blank indicating uncertain), 

marital status of programmer. 
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CARD GB (2A4, A3, 5(D10. 5, Al)) 

DUMMY(3) - hollerith characters not used, 

FIN(I) - zone boundaries increasing order along the x 
axis, 

BCD(I) - flag to indicate end of input if blank, comma 
means to continue. 

Repeat in format (6(D10. 5, Al))if more than five boundaries 
along the x axis are needed. 


CARD GC - same as CARD GB except for y axis. 

Repeat in format (6(DIO. 5, Al)) if more than five boundaries 
along the y axis are needed. 


CARD GD - same as Card GB except for z axis. 

Repeat in format (6(D10. 5, Al)) if more than five boundaries 
along the z axis are needed. 


CARD GE (A4, A2, 315) 

BCD1 - hollerith Z0NE, 

BCD2 - dummy, 

integers which specify the zone as being the 
NXZN0th zone in the x direction, NYZN0th zone 
in the y direction, and NZZN0th zone in the z 
direction. 


CARD GF (2A4, A3, 5(D10. 5, Al)) 

DUMMY(3) - hollerith characters not used. 

FIN(I) - block boundaries in increasing order along the 

x axis, 


NXZN0 

NYXN0 

NZZN0 
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BCD(I) - flag to indicate end of input if blank, comma 
means to continue. 

Repeat in format (6(D10. 5, Al)) if more than five boundaries 
along the x axis are needed. 


CARD GG - same as Card GF except for y axis. 

Repeat in format (6(D10. 5, Al)) if more than five boundaries 
along the y axis are needed. 


CARD GH - same as Card GF except for z axis. 

Repeat in 7 ormat (6(DIO. 5, Al)) if more than five boundaries 
along the z axis are needed. 


CARDS GI to G0 describe the geometry for a block and must be included 
for each block in me zone. 


CARD GI (A4, A2, 315) 


BCD1 - hollerith BL0C, 


BCD2 - dummy, 


NXBND 

NYBND 

NZBND 


integers which specify the block as being the 
NXBNDth in the x direction, the NYBND^h in the 
y direction, and the NZBND in the z direction. 


CARD GJ (3A4, 10(15, Al)) 


NAM2 

DUM(2) 

INP(I) 

BCD(I) 


hollerith MEDI, 
dummy, 

a list of media sector by sector in the block, 

flag to indicate end of input if blank, a comma 
means to continue. 


Continuation with 12(15, Al). 
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CARD GK (3A4, 10(15, 

Al)) 

NAM2 

hollerith SURF, 

DUM(2) 

dummy, 

inp(i) 

a list of quadratic surfaces appearing in the 
block Numbers must appear in the order the 
surfaces are described on Card GQ, 

9CD(I) 

flag to indicate end of input if blank, a comma 
means to continue. 

Continuation of Card GK in 12(15, Al) format is permissible. 

CARD GL (A4, A2, 1813) 

SI 

hollerith SECT, 

DUM 

dummy, 

IND(I) 

the designation of each sector winch describes 
the position of the sector relative to quadratic 
surfaces, 


+1: sector is on positive side of surface, 

-1: sector is on negative side of surface, 


0: surface is not needed to define sector. 

There must be a Card GL for each sector and references to 
quadratic surfaces must be in same order as they are listed 
on Card GQ. 

CARD GM (3A4, 10(15, 

Al)) (omit if NSTAT on Card GA = 2) 

NAM2 

hollerith REGI, 

DUM(2) 

dummy, 


INPflO - a list of regions sector by sector in the block, 

BCD(I) - flag to indicate end of input if blank, a comma 

means to continue. 

Continuation with 12(15, Al) format is permissible. 
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CARD GN (3A4, 10(15, Al)) (omit if NS TAT on Card GA = 2) , 

I 

NAM2 - hollerith SURF, 

DUM(2) . - dummy, 

INP(I) _ same as for Card GK except for region input in- 

BCD(I) " stead of material inpu^. 

CARD GO (A4, A2, 1813) (omit if NSTAT on Card GA = 2) 

SI - hollerith SECT, 

DUM - dummy, 1 

IND(I) - same as for Card GL except for region input in¬ 

stead of material input. 

Repeat\Cards GI to G0 for each block. 

CARD GP (15, 16A4, A2) • 

N0BD - total number of quadratic surfaces in ihe entire 
system, '■ 

DUM (I) - hollerith characters ignored by the code. (Help- 

I ful in identifying input at a later time.) \ 

CARD GQ (4 (DIO. 5\ A4, IX, Al)) 

C0F(J) - coefficient of the term, 

BCD1(J) - hollerith indicating which term of the equation. 

• XJQ, YSQ, ZSQ, XZ, YX, YZ, XY, ZX, YZ, 

X, Y, Z, or blank are the possibilities,. 

BCD2(J) - a flag which indicates the quadratic equation 
continues. Any non-blank character ends the 
fjeld. The next function must start on row card. 

Repeat Caras GQ until all surfaces have been described. 

A sample of the input is shown in Table BUI. 
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05R General Geometry Input Example 














































B. 3 


Combinatorial Geometry Input Instructions 


The combinatorial geometry input data is read by the J0MIN 
subroutine, except for the region volumes VN0R(I), which are read by 
the GTVLIN subroutine whenever IV0PT = 3. For clarity of terminology, 
the terms ’’regions” and "media" have essentially the same meaning as 
in the 05R Geometry Package, but are constructed in a different manner. 
The term ’’zone” is the same as the ’’region” as defined in the original 
combinatorial geometry package. The term "body” has the same mean¬ 
ing as in the original combinatorial geometry package. 


CARD CGA (215, 10X, 10A6) 


IV0PT - option which defines the method by which region 
volumes are determined; if 

IV0PT = 0, volumes set equal to 1., 

IV0PT = 1, concentric sphere volumes are cal¬ 
culated 

* 

IV0PT = 2, slab volumes (1-dim.) are calculated, 
IV0PT = 3, volumes are input by card, 

IDBG - if IDBG > 0, subroutine PR is called to print re¬ 
sults of combinatorial geometry calculations 
during execution. Use only for debugging, 

JTY - alphanumeric title for geometry input (columns 

21-80). 


CARDS CGB (2X, A3, IX, 14, 6E10. 3) 

One set of CGB cards is required for each body and for the END 
card (see Table B IV). Leave columns 1-6 blank on all continu¬ 
ation cards. 


Not operational. 
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ITYPE - specifies body type or END to terminate reading 
of body data (for example B0X, RPP, ARB, etc.). 

• Leave blank for continuation cards, 

IALP - body number assigned by user (all input body 

numbers must form a sequence set beginning at 
1). If left blank, numbers are assigned sequen¬ 
tially. Either assign all or none of the numbers. 
Leave blank for continuation cards, 

FPD(I) - real data required for the given body as shown 
in Table B IV. 


CARDS CGC (2X, A3, 15, 9(A2, 15)) 

Input zone specification cards. One set of cards required for 

each input zone, with input zone numbers being assigned sequen¬ 
tially. 

IALP - IALP must be a nonblank for the first card of 
each set of cards defining an input zone. If 
IALP is blank, this card is treated as a con¬ 
tinuation of the previous zone card, 

IALP = END denotes the end of zone description. 

NAZ - total number of zones that can be entered upon 

leaving any of the bodies defined for this input 
region (some zones may be counted more than 
once). Leave blank for continuation cards for a 
given zone. (If NAZ £ Oon the first card of the 
zone card set, then it is set to 5). This is used 
to allocate blank common, 

Alternate IIBIAS(I) and JTY(I) for all bodies de¬ 
fining this input zone. 

nBIASfl) - specify the "0R M operator if required for the 
JTY(I) body, 

JTY(I) - body number with the (+) or (-) sign as required 
for the zone description 
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CARDS CGD (1415) 


MRIZ(I) - MRIZ(I) is the region number in which the "I** 1 ” 
input zone is contained (I = 1, to the number of 
input zones). Region numbers must be sequen¬ 
tially defined from 1. 


CARDS CGE (1415) 

MMIZ(I) - MMIZOO is the medium number in which the 
ttlthtf i n p U t zone is contained (I = 1, to the 
number of input zones). Medium numbers must 
be sequentially defined from 1. 


CARDS CGF (7E10. 5) (omit if IV0PT ^ 3) 

VN0R(I) - volume of the "I th " region (I = 1 to MXREG, the 
number of regions). 

Note: If ENDRUN is used to obtain collision density and track 
length per unit volume estimate of fluence, then a data 
statement in ENDRUN must give a relationship between 
region and media. In this case only one medium may be 
in a region. 
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TABLE BIV. 

Input Required on CGB Cards for Each Body Type 
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B. 4 


M0RSEC—Cross-Section Module Input Instructions 


CARD XA (20A4) 

Title card for cross sections. This title is also written on tape 
if a processed tape is written; therefore, it is suggested that the 
title be definitive. 


CARD XB (1315) 


Column 

Name 

Description 

5 

NGP 

the number of primary groups for which 
there are cross sections to be stored. 

Should be same as NMGP input in M0RSE, 

10 

NDS 

number of primary downscatters for NGP 
(usually NGP), 

15 

NGG 

number of secondary groups for which there 
are cross sections to be stored, 

20 

NDSG 

number of secondary downscatters for NGG 
(usually NGG), 

25 

INGP 

total number of groups for which cross 
sections are to be input, 

30 

ITBL 

table length, i. e., the number of cross 
sections for each group (usually equal to 
number of downscatters + number of up- 
scatters + 3), 

35 

ISGG 

location of within-group scattering cross 
sections (usually equal to number of up- 
scatters + 4), 

40 

NMED 

number of media for which cross sections 
are to be stored--should be same as MEDIA 
input in M0RSE, 

45 

NELEM 

number of elements for which cross sections 
are to be read, 
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CARD XB (Cont’d.) 


Column 

Name 

Description 

50 

NMEX 

number of mixing operations (elements 
times density operations) to be performed 
(must be ^ 1), 

55 

NC0EF 

number of coefficients for each element, 
including P Q , 

60 

NSCT 

number of discrete angles (usually 

NC0EF/2i ntegra i), 

65 

ISTAT 

flag to store Legendre coefficients if 
greater than zero. 

CARD XC (1115) 


Column 

Name 

Description 

5 

IRDSG + 

switch to print the cross sections as they 
are read if > 0, if < 0 card sequence is not 
checked, 

10 

ISTR + 

- switch to print cross sections as they are 

stored if > 0, 

15 

IFMU f 

switch to print intermediate results of p’s 
calculation if > 0, 

20 

IM0M f ' 

switch to print moments of angular dis¬ 
tribution if > 0, 

25 

IPRIN f 

switch to print angles and probabilities if 
>0, 

30 

IPUN + 

switch to print results of bad Legendre co¬ 
efficients if > 0, 

35 

IDTF f 

switch to signal that input format is DTF-r\ 
format if > 0; otherwise, ANISN format is 
assumed, 


Switches are ignored if DCTAPE < 0. 
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Column Name 


Description 


40 IXTAPE - logical tape unit if binary cross-section 
tape, set equal to 0 if cross sections are 
from cards. If negative, then the pro¬ 
cessed cross sections and clher necessary 
data from a previous run will be read; in 
this case (IXTAPE < 0) no cross sections 
from cards and no mixing cards may be in¬ 
put. The absolute value of IXTAPE is the 
logical tape unit, 

45 JXTAPE - logical tape unit of a processed cross- 

section tape to be written. This processed 
tape will contain the title card, the vari¬ 
ables from common L0CSIG and the per¬ 
tinent cross sections from blank common, 

50 I06RT - logical tape unit of a point cross-section 

tape in 06R format, 

55 IGQPT - last group (M0RSE multigroup structure) 

for which the 06R point cross sections are 
to be used NMGP). 


CARD XD (1415) (omit if IXTAPE * 0) 

Element identifiers for cross-section tape. If element identifiers 
are in same order as elements on tape, the efficiency of the code 
is increased due to fewer tape rewinds. 


CARD XE (omit if IXTAPE i 0) 

Cross sections in ANISN format if IDTF * 0, otherwise, DTF-IV 
format. Cross sections for INGP groups with a table length ITBL 
for NELEM elements each with NC0EF coefficients. 








CARDS XF (215, E10. 5) (NMX cards are required, omit if 
' IXTAPE <0) 

KM - medium number (media numbers from 1 to 

MEDIA, see Card B, must appear on some XF 
card), 

KE - element number occurring in medium KM (neg¬ 

ative value indicates last mixing operation for 
that medium and at least one negative value is 
required for each medium), 

RH0 - density of element KE in medium KM in units of 

atoms/(barn cm). 


CARDS XG (15) (omit if I06RT * 0) 

NXPM - number of point cross-section sets per medium 
found on an 06R tape, 

= 1, total cross section only, 

= 2, total and scattering cross section, 

= 3, total, scattering, and v *fission cross 
section. 


XCHEKR Card (415) Cross sections and cross-section input data may 
be checked independently of M0RSE utilizing XCHEKR. ( 5 ) The 
input to XCHEKR consists of the cross-section cards XA through 
XG preceded by a card as follows: 


IADJM 

MEDIA 

NMGP 

NMTG 


set greater than zero for an adjoint problem, 

number of cross-section media; should equal 
NMED on Card XB, 

number of primary particle energy groups for 
which cross sections are to be stored; should 
equal NGP (or NGG if NGP = 0) on Card XB, 

total number of energy groups for which cross 
section are to be stored. Should be equal to 
NGP + NGG on Card XB. 
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B. 5 SAMB0 Analysis Input Ins Auctions 


The following data are read from cards by SC0RIN: 


CARD AA (20A4) 

Title information—will be immediately output. 


CARD AB (815) 


ND 

- number of detectors (set =1 if £ 0), 

NNE 

number of primary particle energy bins to be 
used, (must be < NE), 

NE 

- total number of energy bins (set =0 if £ 1), 

NT 

- number of time bins for each detector (may be 
negative, in which case | NT j values are to be 
read and used for every detector) (set =0 if 

1nt| s .i). 

NA 

number of angle bins (set =0 if £ 1), 

NRESP 

- number of energy-dependent response functions 
to be used (set =1 if £ 0), 

NEX 

number of extra arrays of size NMTG to be set 
aside (useful, for example, as a place to store 
an array of group-to-group transfer probabilities 
for estimator routines. If the subroutine 
ENDRUN which outputs fluence estimates from 
collision and track lengths is used, then this 
number must be at least MXREG + 2 (see Card I 
for MXREG), 

NEXND 

number of extra arrays of size ND to be set 
aside (useful, for example, as a place to store 
detector-dependent counters). 
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CARDS AC (3E10.4) (ND cards will be read) 

X, Y, Z - detector location. (If other than point detectors 
are desired, the point locations must still be in¬ 
put and can be combined with additional data 
built in to user routines to fully define each de¬ 
tector. ) 

Note that the distanc - between the above points and the XSTRT, 
YSTRT, ZSTRT values and the initial age, AGSTRT, will be 
used to define the lower limit of the first time bin. 


CARD AD (20A4) 

Title or units for total responses for all detectors. Will be used 
in columns 54 through 133 of the title for the print of these arrays. 


CARD AE (20A4) 

Title or units for each total response for all detectors. 


CARDS AF (7E10.4) 

Response function values. NMTG values will be read in each set 
of AF cards. Input order is from energy group 1 to NMTG (order 
of decreasing energy). 

Note: Cards AE and AF are read in the following order: AEj, 
AFj,..., AF n , AE2 , AFj, ..., AF n , etc. NRESP sets of 
AE, AF cards will be read. 


CARD AG (20A4) (omit if NE s 1) 

Units for energy-dependent fluence for all detectors. 


157 





CARDS AH (1415) (omit il NE * 1) \ 

Energy group ^umbers defining lower limit of energy bins (in 
order of increasing group numberJT"The NNE (if >0) entry must 
equal NGPQTN; the NE entry must be set to NMGP ■; NGPQTG 
for a combined problem, or else NGPQTG or NGPQTN.' 

CARD AI (20A4) (omit if | NT) s, 1) ' 

I \ 

Units for time-dependent total responses for all detectors. ' 

CARD AJ (20A4) (omit if | NT j £ 1 or jNE £ 1) ' 

Units for time and energy-dependent fluence for all detectors 

I 

CARDS AK (7E10. 4) (omit if | NT | s 1) ' 

NT vai-es of upper limits of time bins for each detector (in order 
Lit detector number). ■ The values for each 

detector must start on a new card. | NT f values only are read 
if NT is negative. They are then used for every detector. 

CARD AL (20A4) (omit if NA s 1 ) 

! 

Units for angle- and energy-dependent fluence for all detectors. 

CARD AM (7E10.4) (omit if NA * 1) 

Si IS** 5 upper limits o{ bins (actually cosine bins; 

the NA 1 ** value must equal one). 

Following the input for the SAMB0 analysis module, input cards 
for user-written routines INSC0R, SOURCE, and ENDKUN. 




6 Tnput to PICTURE 


Card PA: 
NUSE: 


Format (15) , 

The number of characters to reald on Card B to replace 
the standard values of ATABLE. Leave Card A blank 
and omit Card B if the standard ATABLE is desired. 
NUSE * 50. 


Card PB: Format (50A1) (omit if NUSE = 0) 


ATABLE (I), I = 1, NUSE: The list of characters that are to be 
printed for each medium. For medium N, ATABLE 
(N+l) is printed. \ If N £ 47, ATABLE (48) is printed. 
The standard values of ATABLE are: 

I 

Medium Number Character Printed 

0 (external void) 

1 through 9 1 through 9 

10 through 35 | A through Z 


36 through 46 various special. 

characters 


\ ^ 47 (including 

I internal voids) (blank) 


GE0M input: Combinatorial geometry input. 

C£rd PC: Format (212, 18A4) 

= 0 After this picture, return to Card PC for another 
picture with the same geometry. 

= 1 After this picture, return to Card PA to read in a 
ne’" GE0M input. 

= -1 Display the region geometry. 

IRG: <? = 0 Display the zone geometry. 

( = 1 Display the material geometry. 

'(TITLE(I), ^1 = 1, 18: 72 characters to be printed as a title. 


ICNT: 
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X LR ) 

Y LR / X ’ Y > Z coorc ^ na ^ es in the combinatorial geometry 
( of the lower right corner of the picture. 

Z LR ) 


NOTE: Card PD partially describes the plane of the slice by 

defining two points in the plane and designates the top, 
bottom, left and right sides of the picture. 


6 . 


Card PE: Format (6E10. 5) 


Direction numbers proportional to the direction cosines 
for the U axis of the picture. The U axis points down the 
printed page in the direction the page moves through the 
printer. 



Direction numbers for the V axis of the picture. The V 
axis points to the right across the page. 


NOTE: Card PE completes the description of the plane of the 

slice by giving a line in the plane, also specifies the 
orientation of the picture on the output. 


7. Card PF: Format (215, 2E10. 5) 

NU: Number of intervals to print along the U axis (overrides 

DELU) . 

NV: Number of intervals to print clong the V axis (overrides 

DELV) . 
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DELU: Spacing x a/i GE0M units) of intervals along the U axis. 
DELV: Spacing (in GE0M units) of intervals along the V axis. 

NOTE: All four entries are not required as input on Card PF, 

see below for explanation. 


Options: 


For this case NU and NV must be specified. In addition, either 
DELU or DELV must be specified. If the other is left blank, the 
code will produce an undistorted picture. If both DELU and DELV 
are specified the picture is likely to be distorted. Tfc?' standard 
printers give 10 characters to the inch across a line but only 6 
lines per inch down the page. Because of this DELV = . 6*DELU 
is necessary to produce an undistorted picture. 

X LH 0rY LR 0r W°- 

If any one variable on Card PD is specified, the code will calculate 
the others to produce an undistorted picture. 

If both NU and DELU (or both NV and DELV) are specified, DELU 
(or DELV) will be ignored. 

The U and V axes may have arbitrary orientation. (If they are not 
orthogonal, the resulting picture will be distorted.) In Option 1, 
the first point will be at (X, Y, Z)ul, and the remaining points in 
the directions and at the distances specified. In Option 2, the 
range from Xjjl to X£, R is ^-'vided into intervals and the calcu¬ 
lated points will be at the midpoints of the intervals. The first 
point will be 1/2 interval past (X, Y, Z)ul and final point will 
be within 1/2 interval of (X, Y, Z) LR . If (X, Y, Z) LR does not 
lie on the U-V plane, or if the U and V axes are not orthogonal, 
the location of the final point is not readily predictable. 

The simplest method to obtain the correct results is to specify two 
diagonal corners of the plane of the slice on Card PD, with the top 
having the short dimension, if it is not square, and the side having 
the long dimension. Then, on Card PE, specify the U axis to be 
parallel to the edge of the slice with the large dimension (left or 
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right side), and the V axis to be parallel to the edge of the slice 
with the small dimension (top or bottom). Finally, let the only 
entry on Card PF be NV equal to maximum number of characters 
per line on your printer, this will provide the largest undistorted 
picture. 
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